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by 

Donald  E.  Barrick 


ABSTRACT 


The  effect  of  sea  state  on  ground-wave  propagation  across  the  ocean  is 
computed  in  the  HF  and  VHF  region.  The  history  and  present  understanding  of  ground- 
wave  propagation  is  briefly  reviewed,  especially  as  concerns  the  influence  of 
roughness.  The  approach  of  the  analysis  here  is  to  derive  an  effective  surface 
impedance  at  grazing  which  includes  the  effects  of  roughness.  To  do  this,  the  statis¬ 
tical  boundary  perturbation  approach  of  Rice  is  applied  to  the  sea  surface,  which  is 
"slightly  rough"  at  HF/VHF.  In  addition,  the  Leontovich  (or  impedance)  boundary 
condition  is  employed  because  ocean  water  is  a  good  (but  not  perfect)  conductor 
at  these  frequencies.  The  analysis  shows  that  the  total  effective  impedance  at  the 
surface  can  be  expressed  as  two  terns:  (i)  the  impedance  of  a  perfectly  smooth  sea 
water  surface  at  grazing,  and  (ii)  a  second  term  accounting  for  roughness.  The  latter 
is  obtained  from  the  ocean  wave-height  spectrum. 

The  report  examines  two  height-spectrum  models  for  wind-driven  ocean  waves: 
a  directional  Neumann-Pierson  model  and  an  isotropic  Phillips  spectrum.  The  effective 
surface  impedance  is  calculated  for  these  models.  This  impedance  is  then  used  to 
compute  the  ground-wave  transmission  loss  across  the  sea.  Graphs  are  shown  for  a 
variety  of  frequencies,  ranges,  sea  states,  and  receiver  heights.  Examples,  in  which 
these  curves  are  used  in  communications  problems  are  solved. 

A  bibliography  of  ground-wave  (open  literature)  publications  is  included 
as  an  Appendix.  The  articles  are  arranged  therein  chronologically,  by  decades,  and 
then  alphebetically,  by  author. 
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I .  INTRODUCTION 


Interest  in  ground-wave  propagation  began  with  Sommerfeld' s^  analysis  in 

[2-4] 

1909.  Widespread  utilization  of  these  results,  however,  came  with  Norton's  papers 
which  reduced  Sommerfeld's  complex  expressions  to  graphs  suitable  for  engineering 
applications  and  extended  the  results  to  spherical  surfaces^.  Since  then,  many 
investigators  have  made  contributions  which  include  a  variety  of  deterministic  surface 
geometries,  including  layers,  surface-step  discontinuities,  propagation  beyond  hills 
(of  given  geometrical  shape),  knife-edge  discontinuities,  a  layered  atmosphere,  and 
abrupt  changes  in  surface  properties  (i.e.,  a  shoreline  model).  For  several  recent 
reviews  of  the  pertinent  literature,  one  should  consult  the  articles  by  Wait^, 
Feinberg^,  Bremmer^,  Goubau^  ,  King^1^,  and  Wait^^.* 


With  few  exceptions,  all  of  the  articles  on  ground-wave  propagation  have 

treated  the  earth  surface  as  smooth,  either  planar  or  spherical.  Techniques  have 

been  developed  for  treating  ground-wave  diffraction  by  deterministic  shapes  in  given 

locations,  such  as  a  parabolic-shaped  hill^^,  a  knife-edge ,  a  surface  step1^'*1^, 

[7  15-171 

or  a  linear  interface  between  two  surface  media  (modeling  a  coastline)  ’  .  All 

of  these  obstacles,  however,  are  handled  deterministically,  i.e.,  their  shape  and 

locations  are  specified,  and  the  exact  expression  is  sought  for  their  diffraction 

effect.  None  of  these  surface  irregularities,  either  singly  or  in  combination,  can 

give  a  really  satisfying  model  for  a  statistically  rough  surface,  such  as  the  ocean. 

Hence,  persons  interested  in  ground-wave  propagation  across  the  sea  have  all  used 

Norton's  (or  similar)  techniques  and  assumed  the  sea  surface  to  be  perfectly  smooth. 

* 

Yet  it  has  been  known  for  years  that  a  slight  corrugation  on  an  otherwise  per- 

r  j^g  i 

fectly  conducting  surface  has  a  substantial  effect  on  a  wave  propagating  across  it1 
in  particular,  the  reactive  portion  of  the  surface  impedance  is  enhanced.  Thus  there  is 
good  reason  to  believe  that  in  the  HF  and  VHF  regions,  where  ocean-wave  heights  can 
be  a  significant  percent  of  a  wavelength,  the  roughness  will  play  an  important  role 

and  must  not  be  neglected.  With  a  lack  of  any  theoretical  guidelines,  measurements 

[19] 

carried  out  in  the  past  have  proved  to  be  relatively  meaningless,  especially  since 


*A  novice  to  the  subject  would  do  well  to  consult  one  of  these  reviews,  with  possibly 
a  preliminary  persual  of  a  sound  tutorial  treatment,  such  as  is  contained  in  Jordan's 
textbook! 12}  (Chapter  16).  This  suggestion  is  made  because  hundreds  of  articles  have 
been  published  on  the  subject,  and  any  attempt  to  be  comprehensive  is  bound  to  be 
quite  voluminous.  Rather,  the  reader  may  refer  to  the  Appendix  for  a  chronological 
bibliography  on  most  of  the  important  contributions  on  ground  waves. 
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no  attempts  ware  made  to  ascertain  tha  atatiatlca  of  tha  aurfaca  (i.a.,  eha  aaa  atata) 
at  tha  time  of  tha  measurements.  Sinca  an  ever-increasing  number  of  ayataaa  (both 
communication  and  radar)  will  Involve  propagation  over  the  aaa,  an  effort  to  treat 
tha  effect  of  aurfaca  roughnesa  in  a  quantitative  manner  appear*  in  order. 

There  have  been  a  handful  of  attempt*  to  deal  with  roughnaaa  In  ground-wav* 

propagation  quantitatively.  Feinberg^^  in  an  English  article  showed  already  in 

1944  that  roughnaaa  produces  an  increase  in  the  effective  surface  Impedance,  and 

derived  an  Integral  for  thia  average  impedance.  He  obtained  this  result  from  the 

integral  equation  for  the  fields  obtained  using  Green's  theoram;  in  this  respect, 

Feinberg'a  approach  to  the  problem  is  similar  to  others  who  have  used  integral 

aquations,  e.g.,  Hufford^^  and  King^*^ .  While  Felnberg  initially  applied  the 

Leontovich  boundary  condition,  he  later  presents  the  effective  surface  Impedance  term 

contributed  by  tha  roughness;  the  letter  contains  no  dependence  upon  the  surface 

material  properties,  i.a.,  it  la  strictly  valid  only  in  the  limit  of  a  perfectly 

conducting  surface.  This  useful  relationship,  nonetheless,  apparently  went  unnoticed 

and  unapplied  in  the  West,  and  even  in  the  Soviet  literature  it  was  only  repeated 

(221 

a  couple  of  times*  .  To  tha  best  of  our  knowledge,  no  on*  in  either  ccuntry  has 

F  81 

ever  used  it  to  study  propagation  across  a  rough  sea.  Brenner1  1  later  notes 
Feinberg'a  raault,  and  repeata  it. 

Rice,  on  the  other  hand,  studied  very  thoroughly  the  problem  of  reflection 

*  f  231 

of  electromagnetic  waves  from  a  slightly  rough  surface1  .  Using  the  boundary 
conditions  along  with  Maxwell's  aquations,  he  solves  the  problem  using  a  classical 
perturbation  approach.  In  his  Section  6,  he  briefly  considers  propagation  over  a 
perfectly  conducting  rough  surface  for  vertical  polarization.  His  result  is  identical 
to  Feinberg'a  if  one  makes  tha  connection  between  his  effective  horizontal  propagation 
constant  and  Feinberg'a  surface  impedance.  While  Rice's  classical  paper  has  been 
quoted  and  applied  often  for  scatter  problems,  this  important  result  for  surface-wave 
propagation  has  also  gone  unnoticed. 

[24] 

Walt  considers  propagation  along  a  perfectly  conducting  surface  having 

a  random  distribution  of  hemispheric  bosses  superposed.  These  bosses,  all  of  the  same 

size,  are  small  compared  to  a  wavelength  but  their  spacing  is  large  compared  to  their 

F  251 

radius;  the  work  is  based  upon  multiple  scattering  technique*  derived  by  Twersky  . 

♦"Slightly  rough"  is  the  accepted  term  for  surfaces  whose  roughness  height  is  less 
than  a  wavelength.  The  sea  at  HF  and  much  of  VHF  falls  Into  this  category. 
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The  interpretation  of  these  results  by  Wait  constitutes  possibly  the  first  attempt  to 
consider  the  effect  of  a  random  roughness  on  the  propagation  of  surface  waves.  It  is 
obvious,  however,  that  quantitative  comparison  of  such  a  model  with  the  sea  is  not 
promising.  Senior^ analyzes  the  influence  of  a  slight  roughness  on  the  surface 
impedance,  but  unfortunately  he  does  not  consider  in  sufficient  detail  the  pathological 
case  here  of  grazing  incidence;  his  results  are  mainly  applicable  at  higher  angles  of 
incidence. 

11.  DEFINITIONS  AND  BACKGROUND 

Before  proceeding  further,  it  is  necessary  to  define  terms  to  be  used  here, 

including  ground  wave  and  surface  wave.  The  definition  employed  here  ia  that  of 
(2-5  121 

Norton  ’  .  The  ground  wave  ia  the  total  field  observed  at  a  point  in  space  due 

Co  a  radiating  source  a  finite  distance  away,  excluding  any  component  reflected  from 
the  ionosphere  or  other  discontinuities  in  the  upper  atmosphere;  these  latter 
components  are  termed  sky  waves  and  will  not  be  treated  here.  The  ground  wave  ia 
then  broken  down  further  into  a  space  wave  and  a  surface  wave.  The  space  wave,  if  it 
exists,  consists  of  the  direct  ray  and  that  reflected  from  the  earth;  these  are 
predicted  by  standard  ray  optics  and  exist  only  when  the  observation  point  ia  above 
the  horizon.  The  surface  wave  Is  then  the  remaining  field  in  the  ground  wave  after 
the  space  wave  has  been  subtracted  off.  Persons  using  this  definition  heve  alternately 
referred  to  the  latter  as  the  "Norton  surface  wave". 

It  is  useful  at  this  point  to  review  very  briefly  the  controversy  over  the 
definition  of  surface  waves.  While  we  realize  that  the  reader  of  this  report  is 
little  Interested  in  a  blow-by-blow  description  of  the  quarrel,  it  is  nonetheless 
necessary  because  the  analysis  that  follows  employs  soma  of  the  concepts  subscribed 
to  by  both  camps.  The  controversy  involves  more  than  a  mere  definition  or  terminology; 
it  is  a  question  of  the  nature  of  the  fields  at  the  surface  themselves. 

Sonnet  f eld 's  solution  is  based  upon  a  representation  of  the  radiation  from 
a  dipole  by  a  summation  (l.e.,  an  integral)  of  plane  waves.  This  integral,  along  with 
similar  representations  for  the  reflected  and  transmitted  fields,  satisfy  Maxwell's 
equations  and  the  boundary  conditions  at  the  surface.  The  resulting  expression  for  the 
total  field  above  the  surface  is  then  expanded  asymptotically  (valid  where  distance  from 
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th«  source  la  large  In  terms  of  wavelength).  When  the  terma  appearing  in  this  asymp¬ 
totic  representation  are  interpreted  physically,  the  "space  wave"  tens  described  above 
stand  out;  hence,  it  was  natural  to  call  the  remainder  of  the  solution  a  "surface 
wave".  As  one  approaches  grazing  propagation  on  a  smooth  planar  surface  of  finite 
conductivity,  the  space  wave  vanishes  as  it  should  due  to  the  cancellation  of  the 
direct  and  reflected  waves;  thus  any  power  received  is  due  to  propagation  via  this 
"surface  wave".  Here,  however,  one  can  encounter  difficulties;  one  at  first  night  ex¬ 
pect  that  this  "surface  wave"  alone  should  saLlsfy  Maxwell’s  equations  and  the  boundary 
conditions;  it  does  not.  This  fact  has  caused  doubt  in  the  minds  of  soma,  who  felt  that 
the  entire  expansion  was  somehow  erroneous  or  inadequate.  It  should  have  been  no 
particular  cause  for  alarm,  however.  Many  popular  asymptotic  expansions  of  acattered 
and  diffracted  fields  fall  these  requirements  when  individual  terms  ara  interpreted 
seperetely.  Physical  optics  falls  these  conditions  in  certain  regions  of  space,  and 

yet  it  proves  to  be  a  valuable  tool  when  used  properly.  Individual  creeoing  wave 

I  271 

terma  in  Keller's  geometrical  theory  of  diffraction  fall  these  requirements,  yet 
this  asymptotic  technique  is  considered  an  important  and  meaningful  approximation. 

The  point  to  be  made  here  is  chat  it  is  not  unusual  that  separate  terms  of  a  high- 
frequency  asymptotic  expansion  will  fall  to  satisfy  Maxwell's  equations  and  the 
boundary  conditions,  even  hough  certain  of  these  terms  appear  to  be  separate  physical 
entitles  capable  of  standing  alone.  When  one  employs  all  of  the  terms  in  the  asymptotic 
expansion  together,  and  approaches  the  asymptotic  limit  (usually  allowing  frequency 
to  become  infinite),  the  entire  result  will  nearly  always  satisfy  these  requirements. 

1 28]  * 
Zenneck  first  called  attention  to  the  fact  that  when  one  merely  solves 

Maxwell’s  equations  at  a  planar  interface  between  two  homogeneous  media  (one  of  which 

may  be  slightly  dissipative,  representing  the  ground),  and  applies  the  boundary 

conditions,  one  arrives  at  a  field  which  appears  to  be  "attached"  to  the  surface 

(l.e.,  Its  amplitude  attenuates  exponentially  with  height)  and  which  falls  off 

exponentially  along  the  propagation  direction.  This  solution  satisfies  the  needed 

requirements,  but  assumes  the  media  are  source-free.  The  exponential  attenuation  in 

the  propagation  direction  is  not  entirely  unfamiliar,  and  explains  the  removal  of 

energy  from  the  field;  in  this  sense  it  compiles  with  conservation  of  energy.  This 

removal  of  energy  occurs  here  due  to  the  dissipation  of  the  ground,  which  converts 

electromagnetic  energy  to  heat,  in  other  cases,  such  as  propagation  across  a 

*In  this  re3pect  the  approach  is  similar  to  "assuming"  incident  plane  wav*  fields 
everywhere,  which  is  ccrnmtnly  done  in  certain  situations. 
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slightly  rough  surface  or  propagation  through  rain,  this  attenuation  is  due  to 

removal  and  scatter  of  constant  fractional  increments  of  energy  per  unit  length 

along  the  propagation  direction.  Since  this  wave  obviously  satisfied  the  proper 

equations  to  justify  its  existence  at  an  interface,  Zenneck  felt  Chat  ic  should 

be  called  a  "surface  wave".  Unfortunately,  an  error  in  Sommerfeld's  original  analysis 

resulted  in  a  form  for  Sommerfeld's  surface  wave  which  coincided  with  Zenneck 'a 

definition.  Zenneck  took  this  as  convincing  proof  for  his  explanation  of  radio 

wave  propagation  over  a  surface.  Sommerfeld  corrected  this  oversight  in  a  1926 
[291 

work  ,  which  was  apparently  not  noticed  by  everyone.  By  1936,  serious  questions 

began  to  arise  concerning  the  importance  of  Zenneck's  wave  with  radiation  from 

typical,  finite- length,  radio  antennas;  this  attention  was  spurred  by  experimental 

evidence  which  showed  propagation  behavior  unpredicted  by  the  Zenneck  wave. 

[21  [ 301  [321 

Several  investigations  at  that  time  by  Norton  ,  Rice  ,  and  Burrows  found 

the  error  in  Sommerfeld's  original  work  and  obtained  asymptotic  expansions  for 

the  radiated  field  exhibiting  a  common  behavior.  All  of  these  works,  along  with 

Sonnerfeld 's  1926  paper,  showed  that  if  the  Zenneck  wave  was  present,  it  certainly 

[311 

would  not  be  dominant  at  large  distances  from  the  source.  Yet  they  did  show 
that  the  ratio  between  the  horizontal  and  vertical  components  of  the  electric 
field  at  the  surface  (i.e.,  the  wave  tilt)  is  nearly  identical  to  that  for  the 
Zenneck  wave.  This  ratio  is  proportional  to  Che  normalized  surface, impedance  near 
grazing  Incidence,  and  hence  the  definition  and  use  of  surface  impedance  for  the 

[341 

description  of  ground-wave  propagation  is  valid.  Wait  confirms  this  viewpoint 
and  gives  a  summary  of  ground-wave  propagation  from  a  vertical  dipole  based  upon 

[4] 

the  properties  of  the  surface  Impedance,  Furthermore,  Norton  showed  that  the 
Poyntlng  vector  direction  for  his  surface-wave  component  coincides  nearly  perfectly 
with  that  predicted  for  the  Zenneck  wave  when  the  surface  impedance  is  relatively 
low.  This  Poyntlng  vector,  being  directed  slightly  into  the  surface,  indicates 
that  the  energy  is  flowing  into  the  lower  medium  at  an  angle  very  close  to  the 
Brewster  (or  "pseudo-Brewster")  angle. 
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Confusion  and  controversy  continued.  Several  other  treatments  of  radiation 

from  a  finite  source  were  developed  based  upon  Integral  equations,  a  different  approach 

from  that  used  by  Somnerfeld  and  Norton.  Feinberg^^  and  Hufford^^  employed  the 

Integral  equation  resulting  from  the  application  of  Green's  theorem  to  the  Helmholtz 

[101 

equation  (the  latter  arising  from  Maxwell's  equations).  King  used  a  different 
integral  equation  originating  from  the  "Compensation  Theorem”.  All  of  these  analyses 
are  based  upon  radiation  from  a  finite  source,  e.g.,  a  dipole,  above  a  planar  earth; 
such  a  source  in  the  absence  of  the  plane  would  radiate  a  spherical  wave.  They  all 
arrive  at  Che  same  asymptotic  equation  as  the  Norton-Sommerfeld  result.  Over 
the  years,  experimental  evidence  for  radio-wave  propagation  has  tended  to  confirm 
these  asymptotic  expressions,  rather  than  the  rapidly  attenuating  Zenneck-wave 
explanation.  The  charge  that  the  asymptotic  expansion  for  the  Norton  surface  wave 
does  not  satisfy  Maxwell's  equations  and  the  boundary  conditions  is  answered  in  a 
recent  paper  by  King^*^ .  He  shows  that  the  required  correction  la  arbitrarily  small 
if  the  source  and  observer  points  are  close  to  the  surface  and  the  surface  Impedance 
is  quite  low  (the  latter  condition  implying  an  electrically  dense  and/or  lossy  lower 
medium) .  These  conditions  are  nearly  always  satisfied  in  the  case  of  propagation 
over  the  earth. 

Based  on  the  concurrence  of  the  several  analyses  and  experimental  evidence, 
one  must  conclude  that  the  Norton-Sommerfeld  result  properly  explains  propagation  over 
a  plane  earth  when  the  field  radiates  from  a  finite-sized  source  (e.g.,  a  dipole). 

The  Zenneck  wave  is  not  postulated  on  radiation  from  a  finite-sized  source,  and  hence 

*The  need  for  some  sort  of  "meeting  of  the  minds"  readied  a  culmination  In  the  late 
1950's,  when  a  working  group  was  set  up  by  U.R.S.I.  (International  Scientific  Radio 
Union)  under  the  chairmanship  of  J.  R.  Walt  to  agree  upon  suitable  definitions  on  the 
nature  of  surface  waves.  This  group  more  or  less  dissolved  at  the  URSI  13th  General 
Assembly  in  London,  England,  in  1960  when  it  appeared  that  further  discussion  on  the 
topic  at  that  time  would  not  be  fruitful.  The  Proceedings  of  that  meeting  contain 
discussion  by  participants  reflecting  their  views  on  Che  subject  (pp  533-539);  more 
data  on  the  subject  by  many  of  the  principal  participants  can  be  found  in  the  Trans¬ 
actions  of  the  IRE  on  Antennas  and  Propagation,  Vol  AP-7,  Special  Supplement,  pp  132- 
243A  (1959).  Waiti^M  later  provided  a  needed  comprehensive  historical  review  of  the 
subject;  this  is  both  objective  and  quite  thorough. 

tThe  nature  of  the  surface-wave  over  a  spherical  earth  is  fraught  with  the  same 
general  problems,  i.e.,  does  a  wave  of  the  Zenneck  type  dominate?  Only  the  plane- 
earth  case  is  reviewed  here  for  brevity. 
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1C  does  not  apply  in  the  ground-wave  propagation  problem  since  the  source-generated 
asymptotic  expansions  show  no  evidence  of  its  existence.  Nobody  has  yet  designed  sn 
antenna  which  can  excite  a  Zenneck  wave  over  a  plane  earth  (if  such  a  finite-sized 

i  igi 

antenna  is  even  possible),  although  Barlow  and  Brown1  1  discuss  this  subject. 
Therefore,  one  must  employ  the  source-postulated  approach,  along  with  the  asymptotic 
expansions  attributed  to  Norton  and  Sommerfeld  (both  for  the  plane  and  spherical 
earth)  when  analyzing  ground-wave  propagation  over  the  earth. 

Nonetheless,  paradoxes  remain.  If  the  Zenneck  wave  is  not  an  important 
component  in  the  groundwave  from  a  source,  why  does  the  polarization  or  wave  tilt 
obtained  from  the  asymptotic  solution  agree  so  closely  with  that  of  the  Zenneck 
surface-wave?  Why  does  the  energy  at  the  surface  propagate  into  the  ground  at  the 
same  angle,  i.e.,  the  Brewster  angle,  for  both  cases?  Are  these  facts  merely  a 

[34] 

coincidence?  Wait  points  out  that  rather  than  coincidence,  this  near-equivalence 

of  the  wave  tilt  at  the  surface  is  a  consequence  of  the  fact  that  the  Zenneck-wave 

pole  in  the  Somnerfeld  integral  is  near  the  free-space  wavenumber  In  the  complex 

wavenumber  plane.  This  mathematical  connection  between  the  two,  however,  does  not 

provide  a  clear  physical  interpretation  of  their  relationship.  An  exact  solution  of 

Maxwell's  equations  and  the  boundary  conditions,  in  the  absence  of  sources,  shows  that 

such  a  solution  must  be  a  Zenneck  wave,  which  is  nearly  planar  and  attenuating  in 

both  the  height  and  propagation  directions,  over  a  restricted  localized  region  near 
1 12  18] 

the  surface'  ’  1 ,  Why  is  such  behavior  not  present  at  great  distances  from  the 
source  in  the  asymptotic  results  of  Norton?  Needless  to  say,  until  these  points  are 
answered  in  physically  meaningful  manner,  controversy  and  confusion  is  bound  to 
continue. 


/ 


■7  * 
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111.  APPROACH  AND  SCOPE  OF  PRESENT  ANALYSIS 

The  purpose  of  the  present  analysis  is  to  quantitatively  establish  the 
effect  of  sea  state  on  ground-wave  propagation  and  attenuation  over  the  ocean.  Due 
to  the  nature  of  the  approximations  made,  the  results  will  be  valid  at  frequencies  in 
and  below  the  VllF  region;  vertical  polarization  only  is  considered. 

The  first  portion  of  the  analysis  will  derive  an  expression  for  the  effective 
(or  average)  surface  impedance  of  a  roughened  surface.  Instead  of  treating  the 
surface  as  a  perfect  conductor,  we  shall  use  the  I.eontovich  boundary  condition  for 
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good  (but  finitely)  conducting  surfaces;  this  condition  is  especially  applicable  to 

ocean  water  in  the  HF  and  VHF  regions.  Even  though  the  sea  surface  nay  appear  to  be 

a  good  conductor,  the  distinction  between  a  good  conductor  and  a  perfect  one  at  these 

frequencies  is  an  important  one.  A  ground  wave  propagating  across  a  smooth  perfectly 

conducting  surface  suffers  no  attenuation  due  to  the  presence  of  the  surface,  and 

propagates  as  though  it  were  in  free  space.  A  good  (but  not  perfect)  conducting 

* 

smooth  surface  will  produce  attenuation  of  the  ground  wave  at  the  surface  ;  for  sea 
water,  this  attenuation  factor  results  in  loss  in  signal  level  of  tans  of  decibels  at 
ranges  of  50  miles  or  greater  even  if  the  sea  is  treated  as  perfectly  smooth.  Hence, 
the  present  analysis  of  roughness  will  take  into  account  this  finite  conductivity 
from  Che  outset.  The  Leontovlch  boundary  condition,  in  contrast  with  other  approaches, 
is  discussed  In  the  next  section. 

i  261 

As  shown  in  Senior's  analysis1  ,  the  effective  surface  Impedance  for  any 
surface  is  a  function  both  of  the  angle  of  incidence  and  polarization.  Since  only 
vertical  polarization  is  considered,  the  remaining  question  arises  as  to  the 
angle  of  incidence  to  be  used  here.  This  is  Important,  because  results  useful  for 
scatter  at  higher  Incidence  angles  are  not  usually  meaningful  when  this  angle  is 
reduced  to  the  near-grazing  region  common  in  ground-wave  propagation,  and  one  must 
treat  this  case  separately.  To  do  this,  a  plane  wave  will  be  assumed  to  ba  guided  by 
the  surface,  and  the  effective  vertical  wave  number  will  be  determined  and  related 

to  the  surface  impedanca.  In  this  respect,  the  approach  taken  is  the  seme  as  that 
1231 

in  Rice1  (i.e.,  a  perturbation  approach),  except  for  the  difference  in  boundary 
conditions.  This  guided  wave  is  identical  to  the  Zenneck  wave  discussed  previously. 
One  may  ask  why  such  an  approach  has  validity  here  when  we  ultimately  want  to 
consider  a  finite  source  which  produces  no  Zenneck  wave  in  the  radiated  far  field? 

The  Justification  ia  based  upon  the  discussion  in  the  preceding  section;  there  it 
was  mentioned  that  results  of  Norton^  and  Wise^^  show  that  the  Norton  surface-wave 
derived  from  the  finite  source  radiator,  has  the  same  phase  and  polarization  tilt 
(i.e.,  direction  of  Incidence)  as  the  Zenneck  wave,  even  though  the  attenuation  in 
no  way  resembles  that  of  the  Zenneck  wave.  Hence,  the  affective  surface  Impedance 

derived  by  employing  a  Zenneck  guided  wave,  being  applicable  for  the  same  incidence 

121  1  33  341 

angle,  can  be  used  in  the  source-based  formulation  of  Norton  and  Wait  ’  .  The 

*Thia  has  been  defined  in  the  preceding  section  as  the  Norton  surface  wave. 
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reasons  chat  we  Initially  assume  a  plane  wave  propagating  in  a  guided  mode  (Zenneck 

wave)  across  the  surface  are  primarily  for  the  mathematical  convenience  and  the 

physical  insight  obtainable.  A  check  on  the  validity  is  obtained  by  comparison  of 

(221 

result  with  that  of  Feinberg*  where  we  permit  the  conductivity  to  approach 
infinity.  Feinberg' s  result  was  derived  from  an  integral  equation  approach  involving 
the  source.  Besides  providing  an  alternate  derivation  from  Feinberg' s,  the  analysis 
here  shows  explicitly  the  effect  of  the  finite  conductivity  in  the  roughness  contri¬ 
bution  to  surface  impedance.  This  addition  indicates  the  limits  of  validity  in 
Feinberg' s  assumption,  and  also  yields  a  correction  factor  which  accounts  for  this 
deviation  from  perfect  conductivity.  Furthermore,  it  obviates  a  numerical  difficulty 
encountered  in  computing  the  integral  for  a  given  sea  surface  model. 

With  the  above  technique,  we  derive  an  expression  for  the  average,  or 
effective,  surface  impedance.  This  is  the  sum  of  two  terms:  the  first  being  the 
normal  impedance  of  the  otherwise  smooth  surface,  and  the  second  term  representing 
the  additional  effect  of  the  roughness.  The  scale  and  spatial  height  spectrum  of  the 
roughness  (or  the  roughness  wavelengths)  determine  the  nature  of  this  impedance 
change.  In  the  HF/VHF  region,  the  radio  wavelengths  are  of  the  same  order  as  the 
typical  ocean  wavelengths,  and  a  noticeable  effect  on  the  impedance  is  expected.  In 
order  to  study  this  in  a  more  quantitative  manner,  we  select  two  semi-empirical  models 
for  the  sea  surface:  the  Neumann-Pieraon  spectrum  and  the  Phillips  spectrum.  In  the 
first  model,  the  ocean  waves  are  assumed  to  be  slightly  directional  in  nature,  while 
in  the  latter,  a  non-directlonal  (or  isotropic)  angular  dependence  is  assumed.  These 
models  express  the  ocean  wave-height  strength  in  terms  of  wind  speed,  assuming  that 

this  wind  has  been  blowing  sufficiently  long  that  the  resulting  sea  state  is  fully 
[35] 

developed  ,  Using  the  models,  values  of  effective  surface  impedance  are  computed 
numerically.  This  is  done  for  a  variety  of  frequencies  from  1  to  300  MHz,  and  for 
five  sea  states. 

The  reason  for  computing  the  effective  impedance  of  the  rough  sea  is  to 

eventually  calculate  the  strength  of  the  field  in  the  ground  wave  radiated  from  a 

source  above  such  a  surface.  This  is  done  using  the  computer  program  of  Berry  and 
[36] 

Chrisman  for  ground-wave  propagation  over  a  spherical  earth  for  dipole  antennas. 

In  other  words,  after  the  average  surface  Impedance  is  calculated  whicli  includes  the 
roughness  effects,  we  assume'' that  the  ground  consists  of  a  smooth  surface  with  this 
equivalent  Impedance.  This  procedure  is  in  the  spirit  of  Wait's  article^ ^  which 
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suggest*  this,  for  exsaple,  for  corrugated  surfaces.  Curves  plotted  from  this  pro¬ 
gram  are  Included  which  show  the  basic  transmission  loss  to  various  points  above  the 
earth  as  a  function  of  sea  state,  frequency,  range,  and  height.  This  la  done  for 
propagation  in  both  the  downwind  and  crosswind  directions  so  that  some  variance  of 
signal  strength  can  be  obtained. 

Another  effect  occurs  as  a  result  of  the  roughness.  The  ground  wav*  prop¬ 
agating  across  the  sea  interacts  with  the  ocean  waves,  and  a  portion  of  the  energy  in 
it  Is  removed  and  scattered  off  Into  the  sky.  This  scatter  does  not  take  place  for 
smooth  surfaces.  Hence  the  total  field  arriving  at  any  point  above  the  surface  is 
made  up  of  two  components:  the  ground  wave  which  would  exist  at  this  point  for  a 

perfectly  smooth  surface  with  the  equivalent  impedance  calculated,  along  with  the 

★ 

scatter  from  all  points  on  the  surface  due  to  roughness  .  The  former  component  is 
coherent,  meaning  chat  a  CW  transmitted  signal  will  arrive  at  this  point  as  a  pure 
CM  signal.  The  latter  scatter,  however,  is  incoherent,  in  that  its  spectrum  will  no 
longer  be  a  pure  CW  sine  wave,  but  will  be  modulated  by  the  time-varying  ocean 
surface.  The  spectral  spreading  of  the  latter  can  be  of  the  order  of  several  tenths 
of  a  hertz.  Curves  of  Che  total  average  intensity  of  this  incoherent  component  will 
be  shown  in  a  subsequent  report  and  compared  to  the  coherent  power  In  the  ground  wave 
below  it  at  the  sea  surface. 

Additional  subjects  such  as  the  sea  clutter  power,  its  spectrum,  and  bi- 
static  geometries,  will  be  deferred  to  a  later  report.  A  section  will  describe  the 
use  of  the  transmission  loss  curves  for  ccmmnmications  applications. 


IV.  ANALYSIS 

A.  Formulation  of  Various  Boundary  Conditions 

1.  General  Conditions  at  Interface 
Between  Two  Homogeneous  Media 

To  solve  any  problem  involving  Che  interaction  of  an  electromagnetic  wave 
with  an  interface  between  two  media,  one  must  relate  the  field  quantities  at  the 

*This  does  not  take  into  account  any  lonospher leal iy-ret lected  components. 
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interface  by  means  of  constraints  commonly  referred  to  as  boundary  conditions.  The 
most  basic  and  general  conditions  are  derived  directly  from  Maxwell's  equations  in 
their  integral  form.  These  conditions  are  stated  as  follows:  (i)  the  tangential 
components  of  the  electric  and  magnetic  (E  and  H)  fields  uniat  be  continuous  across 
the  interface,  and  (ii)  the  normal  components  of  the  electric  and  magnetic  flux 
densities  (D  and  B)  must  be  continuous  across  the  interface.  This  boundary  condition 
is  valid  for  any  arbitrarily  curving  Interface.  Unfortunately,  it  involves  explicitly 
the  fields  on  both  sides  of  the  interface  (except  when  the  surface  material  is  a 
perfect  conductor;  this  will  be  discussed  in  the  next  section),  which  increases  the 
complexity  of  solving  many  typical  scatter  problems. 

2.  Perfectly  Conducting  Surface 

When  the  surface  material  is  a  perfect  conductor,  the  above  boundary 
conditions  simplify  to  the  following:  (1)  the  total  tangential  electric  field  above 
Che  surface  is  zero,  or  in  equation  form, 

E1  -  (n  •  i1)/!  -  0  (1) 

and  (il)  the  total  normal  magnetic  field  at  the  surface  is  zero,  or  in  equation  form, 

n  •  H*  -  0  (2) 

a  -T 

In  the  above  equations,  n  is  the  unit  normal  to  the  surface,  E  and  H  are  the  total 
electric  and  magnetic  fields  at  the  surface.  Again,  these  conditions  are  valid  at 
any  perfectly  conducting  surface,  regardless  of  whether  it  is  curved  or  planar,  and 
independent  of  Che  angle  of  incidence.  They  are  more  convenient  usually  than  those 
in  (1)  because  they  Involve  only  the  fields  above  the  surface,  in  the  medium  of  Interest. 
They  do,  however,  involve  only  tangential  E  and  normal  H;  i.e.,  there  are  no  similar, 
simple  relations  for  normal  E  and  tangential  H  at  a  perfect  conductor.  To  obtain  such 
relationships  in  terms  of  the  latter  two  components,  one  must  solve  Maxwell's 
equations,  and  the  results  are  not  especially  tractable  or  useful. 


!l-A8LE  COP)' 
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3.  Planar  Interface  Between  Homogeneous 
Media— Fresnel  Reflection  Coefficients 


One  case  In  which  a  simple  boundary  condition  obtains  for  general  homogeneous 
materials  occurs  where  the  interface  is  perfectly  smooth  and  planar.  In  this  case  the 
total  field  on  one  side  of  the  surface  can  be  expressed  entirely  in  terms  of  the 
Fresnel  reflection  coefficients  and  the  incident  (locally  plane)  wave.  These 
reflection  coefficients  are  functions  of  the  incidence  angle  and  surface  material 
constants.  Any  dissipative  loss  in  the  material  (either  electric  or  magnetic)  can  be 
included  in  these  coefficients  as  a  complex  permittivity  and/or  permeability.  The 
elimination  of  the  fields  on  the  other  side  of  the  surface  in  this  boundary  condition 
is  convenient,  but  unfortunately  the  entire  method  is  strictly  valid  otily  for  perfectly 
smooth,  planar  interfaces  (surfaces).  This  is  in  contrast  to  the  boundary  condition 
for  the  perfectly  conducting  surface  discussed  previously,  which  is  valid  for  surfaces 
with  any  degree  of  curvature.  That  boundary  condition  applied  only  to  the  total 
tangential  electric  field  at  the  surface,  whereas  the  one  discussed  here  for  planar 
surfaces  applies  to  the  total  (both  normal  and  tangential  components)  electric  and 
magnetic  fields  at  the  surface.  How  to  deal  with  curving  surfaces  will  be  considered 
subsequently. 

4.  Curving  Interface-Tangent  Plane 

Approximation  and  Fresnel  Coefficients 

The  preceding  section  showed  that  the  total  fields  at  a  perfectly  planar 
surface  can  be  written  directly  in  terms  of  the  Incident  field  times  an  expression 
involving  the  Fresnel  reflection  coefficients.  The  question  arises:  under  what 
circumstances  can  such  a  convenient  expression  be  extended  to  curving  and  rough 
surfaces?  This  question  has  been  examined  in  many  places,  and  the  general  conclusions 
are  that  it  is  valid  when  applied  to  surfaces  whose  local  radii  of  curvature  are  much 
greater  than  wavelength.  In  fact,  the  correction  term  in  this  approximation  can  be 
shown  to  be  of  the  order  of  (kop)-1,  where  ko  *  (X  is  the  free-spaca  wavelength), 
and  p  is  a  radius  of  curvature  of  the  surface.  When  wavelength  and  surface  curvatures 
are  such  chat  the  approximation  Is  reasonably  valid,  one  writes  the  total  field  at  a 
given  point  on  the  surface  in  terms  of  the  incident  field  and  the  reflection  coeffi¬ 
cients,  the  latter  being  taken  at  the  angle  of  Incidence  to  that  point.  Hence  the 
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resulting  expression  for  the  total  field  at  the  surface  varies  from  point  to  point 
with  the  local  incidence  angle;  this  local  incidence  angle  is  defined  as  the  angle 
between  the  incident  wave  propagation  direction  and  the  normal  to  the  surface  at  that 
point.  The  formulation  of  the  total  fields  at  the  surface  in  such  a  manner  is  called 
the  tangent-plane  approximation;  it  is  often  used  along  with  physical  and  geometrical 
optics  methods  to  compute  the  fields  reflected  and  scattered  from  curving  surfaces 
in  the  high-frequency  limit. 

The  tangent-plane  approximation  carries  two  further  restrictions  with  it. 
First,  at  regions  of  the  surface  not  directly  visible  to  the  incident  wave  (i.e., 
shadowed  regions),  the  fields  are  assumed  to  be  identically  zero.  This  vanishing  of 
the  fields  in  shadowed  regions  is  not  strictly  valid  in  itself.  Furthermore,  it 
introduces  complications  in  having  to  determine  (either  deterministically  or  statis¬ 
tically)  the  regions  shadowed  for  different  incidence  angles.  Near  grazing,  this 
shadowing  effect  becomes  so  serious  that  its  neglect  cannot  be  tolerated;  yet  it  is 
at  near-grazing  angles  that  surface  waves  appear  to  propagate.  A  second  restriction 
inherent  in  the  tangent-plane  approximation  is  the  lack  of  multiple  scattering  between 
two  or  more  surface  points.  This  effect  is  neglected  because  the  total  fields  at  each 
surface  point  are  written  in  terms  of  the  incident  field  only  at  that  point,  which 
excludes  any  field  contribution  reflected  from  a  nearby  surface  point.  Multiple 
scattering  is  expected  to  be  a  serious  contributor  only  for  surfaces  having  many  large 
concave  regions,  and  is  not  typical  of  the  sea  surface,  for  example. 


5.  Impedance  (Leontovlch)  Boundary  Condition-: 
Curving  Surface 


A  further  boundary  condition  which  can  be  very  usef»  i  for  certain  situations 
is  called  an  impedance  boundary  condition,  in  particular,  wh  .  -  applied  to  vectorial 
electromagnetic  waves,  it  is  commonly  referred  to  as  the  U'm  >vich  boundary 
condition  because  of  the  detailed  pioneering  invest igation*  /  Leontovich  on 

this  subject  in  the  1930's.  He  showed  that  when  the  condv  .4vity  and  refractive 
index  of  the  medium  below  Che  interface  are  relatively  large,  one  can  arrive  at  an 
expression  for  the  total  tangential  fields  immediately  above  the  surface  which  does 
not  involve  fields  in  the  medium  at  all.  Furthermore,  they  do  not  require  serious 
restrictions  on  the  surface  radii  of  curvature,  as  does  the  tangent-plane  approximation 
discuased  in  the  preceding  section. 
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In  equation  font,  the  Leontovlch  boundary  condition  appears  to  ba  a  alight 
extension  of  Equation  (1)  for  the  perfectly  conducting  surface: 

i1  -  (n  •  E^n  “  ZB(n  x  H*)  .  (3) 

--T  * 

Here,  E  and  H  are  the  total  electric  and  magnetic  fields  above  the  surface  and  n  Is 

the  unit  normal  to  this  surface.  The  quantity  Z#  is  termed  the  surface  impedance. 

For  a  homogeneous  medium  below  the  surface,  Z(  is  given  by 

*.  •  7T  •  <‘> 


where  end  s  are  the  absolute  permeability  and  permittivity  of  the  material. 


Obviously,  the  above  boundary  condition  carries  certain  restrictions  on  the 
media  and  interfaces  to  which  it  applies.  As  derived  by  Leontovlch,  these 
restrictions  are  the  following:  (i)  the  index  of  refraction  (i.e.,  /£*)  of  the 
material  is  large  and  has  a  large  imaginary  part,  (ii)  the  fields  immediately  above 
the  surface  vary  slowly  along  the  surface  over  a  distance  of  the  order  of  a  wavelength 
in  the  material  (i.e.,  a  wavelength  in  the  material  is  X^  *  where  X„  Is  the 

wavelength  in  free  space),  and  (ill)  the  radii  of  curvature  of  the  surface  are  small 
compared  to  X  ,  the  wavelength  in  the  material. 

fll 

Thus,  while  there  is  a  restriction  on  the  radius  of  curvature  of  the  surface, 

it  is  considerably  less  stringent  than  that  for  the  tangent-plane  approximation.  The 

condition  there  was  that  the  radii  of  curvature  had  to  be  large  compared  to  X<j,  the 

[  381 

wavelength  in  free  space.  Leontovlch1  1  in  fact  derives  the  firq,t-order  correction 
tenns  for  the  above  boundary  condition  (Equation  (3))  and  finds  that  for  a  homogeneous 
material,  they  contribute  terms  to  the  right  side  of  (3)  of  the  order  of 


(3) 


where  pi  and  pa  are  the  principal  radii  of  curvature  of  the  surface  at  any  point. 


Investigators  have  found  that  the  impedance  boundary  condition  expressed 
in  (3)  is  useful  not  only  for  homogeneous  media  below  the  surface,  but  also  for  other 
types  of  surfaces  such  as  corrugated  and  dielectric-clad  conducting  surfaces.  In 
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fact,  the  purpose  o£  this  report  Is  to  derive  an  effective  Z  for  a  randomly  rough 

•  ® 

sea  surface.  A  rather  thorough  review  of  the  Leontovlch  boundary  condition  and  the 

general  impedance  boundary  condition  for  various  types  of  surface  conditions  was 

.  .  .  .  {39,26] 

done  by  Senior  . 

B.  Applicability  of  Boundary  Conditions  to  the  Sea  at  I1F/VHF 

Not  all  of  the  above  boundary  conditions  are  applicable  to  a  rough  sea  in 

the  HF/VHF  region.  The  boundary  condition  selected  depends  upon  the  method  of  attack 

used  to  solve  the  problem  and  upon  the  material  and  surface  properties  of  the  ocean. 

Several  investigators  have  employed  questionable  (and  even  inapplicable)  boundary 

conditions  and  techniques  in  the  past  to  study  sea  scatter  at  HF/VHF,  and  this  is  the 

reason  that  the  question  is  studied  rather  thoroughly  here.  It  is  the  conclusion  of 

this  section  that  the  only  strictly  valid  technique  to  determine  scatter  from,  and 

propagation  along,  a  rough  sea  at  HF/VHF  is  the  boundary  perturbation  approach 
[231 

attributed  to  Rice  .  Using  this  technique,  the  only  strictly  applicable  boundary 
condition  is  the  Leontovlch  condition,  for  scatter  and  propagation  at  near-grazing 
Incidence. 


1.  General  Conditions  at  Interface  Between 
Two  Homogeneous  Media 

Rice  initially  formulated  a  solution  for  scatter  from  a  slightly  rough 

random  surface  between  two  dielectric  media  for  a  horizontally-polarized  incident 

w*ve.  His  results  were  extended  to  vertical  incident  polarization  by  Peake^0^ ,  and 
[41  421 

by  Barrick  and  Peake  *  to  homogeneous  surfaces  with  arbitrary  g  as  well  as  «. 

All  of  these  analyses  are  based  upoi.  a  perturbation  expansion  of  both  thn  surface  and 
the  reflected  fields.  The  boundary  condition  used  is  the  general  condition  involving 
the  fields  on  both  sides  of  the  interface,  as  discussed  under  1.  of  the  preceding 
section.  While  the  expansion  is  messy  algebraically,  it  is  straightforward,  and 
results  may  be  found  in  the  above  references. 

As  applied  to  the  sea  surface  at  HF/VHF,  these  results  are  believed  to  be 
very  questionable,  if  not  entirely  Invalid.  This  is  due  to  the  fact  that  tiie 
expansion  of  the  perturbed  fields  beneath  the  surface  is  based  upon  a  simple  series 
representation  of  the  quantity  in  which  only  first  three  terms  are 
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retained.  This  is  valid  if  c(m,n)C  is  small  (£  is  tha  ocaan  wava  height  above  a 

\t 

mean  level).  However,  0(10,0)  is  defined  as  V 

*  a  • 

c(m,n)  -  ko/*r>*r  "  “  sin^G^  ,  (6) 

where  9#  is  Che  scattering  angle  from  tha  vertical  to  the  earth,  and  c^,  are  the 

relative  constitutive  parameters,  which  in  general  may  be  complex.  For  tha  see  the 
dielectric  constant  is  about  80  and  the  average  conductivity  is  about  4  mhos/m.  At 
10  MHz,  this  gives  >1,  -  80  +  17200.  For  real  scattering  angles  8^,  |c(m,n)| 

is  of  the  order  of  85  ko  *■  18  m"1.  Hence,  even  for  a  very  calm  sea  with  ocean  waves 
whose  heights,  £,  are  of  the  order  of  6in,  the  quantity  |c(m,n)C|  is  of  the  order  of  3, 
and  clearly,  the  first  three  terms  of  Che  exponential  alc(°»nK  are  entirely  insuf¬ 
ficient.  In  addition,  the  representation  of  the  solution  with  a  perturbation  formula¬ 
tion,  where  c(m,n)£  is  a  "smallness"  parameter,  is  invalid  because  such  a  smallness 

* 

parameter  must  always  be  considerably  less  than  unity.  The  above  restriction  on  ocean 
waves  fails  throughout  Che  entire  HF/VHF  region. 

2.  Perfectly  Conducting  Surface 

The  ocean  at  HF/VHF  is  a  good  conductor,  and  for  many  applications  may  be 
considered  a  perfect  conductor.  For  our  application  here,  however,  such  a 
simplification  ia  not  possible.  The  following  is  an  explanation  of  the  reason.  A 
ground  wave  propagates  in  a  direction  close  to  grazing,  but  ac'.ually  appears 
to  be  propagating  into  the  surface  at  the  Brewster  angle.  For  sea  water,  this  angle 
at  10  MHz  is  about  1/2  degree  from  grazing.  For  the  vertical  polarization  states, 
propagation  across  a  perfectly  smooth,  perfectly  conducting  surface  takes  place  with 
no  attenuation  (other  than  that  of  free  space).  In  addition,  scatter  from  a 
perfectly  conducting,  slightly  rough  surface  is  entirely  insensitive  to  grazing 
angle  in  this  region,  i.e.,  the  scattered  pwer  for  incidence  and  scattering  angles 
of  about  2#,  1/2*,  and  0*  above  grazing  are  identical. 

*The  discussion  here  is  not  meant  to  imply  that  Rice's  analysis  using  the  interface 
boundary  conditions  between  two  homoegenous  media  is  invalid  in  general.  There  are 
a  great  many  surfaces  and  materials  which  meet  this  restriction.  The  condition  is 
easier  to  satisfy  for  less  dense  surface  materials;  sea  water,  however,  is  quite  dense 
electromagnetically  at  HF/VHF. 
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For  a  good,  but  not  p«r£ect,  conductor  however,  propagation  near  grazing  in 

the  vertical  polarization  states  differs  considerably  from  the  behavior  described 

above.  First,  ground-wave  propagation  across  a  perfectly  smooth  sea  surface  will 

suffer  several  tens  of  decibels  attenuation  over  that  a  perfectly  conducting  surface 

at  a  range  of,  say,  100  miles.  Hence,  by  assuming  a  perfectly  conducting  sea,  one 

loses  entirely  the  ground-wave  attenuation  function  which  is  the  essence  of  the 

problem.  Second,  when  such  a  surface  is  slightly  rough,  the  scatter  for  the  vertical 

polarization  states  has  been  shown  to  vary  considerably  at  angles  near  grazing  , 

★ 

in  contrast  with  the  case  of  a  perfectly  conducting  surface.  Here  again,  the  dif¬ 
ference  between  a  good  and  a  perfect  conductor  is  critical. 

For  HF/VHF,  and  for  most  sea  scatter  problems  involving  horizontal  and 

vertical  polarization  (where  incidence  and  scatter  angles  are  restricted  from  the 

region  within  3“  of  grazing),  the  assumption  of  a  perfectly  conducting  surface  is 

valid.  In  these  cases,  the  Rice  perturbation  technique  based  upon  equation  (1)  is 

entirely  adequate,  and  results  derived  and  presented  for  average  scattering  cross 
T  40-421 

sections  in  are  applicable  (l.e.,  where  a  perfectly  conducting  surface  is 

assumed).  For  the  near-grazing  angles  associated  with  ground-wave  propagation,  however, 
the  assumption  of  perfect  conductivity  can  lead  to  dangerous  and  possibly  erroneous 
results  and  interpretations,  as  seen  from  the  discussion  in  the  preceding  paragraph. 

3.  The  Tangent-Plane  Approximation 

The  tangent-plane  approximation,  as  discussed  in  the  preceding  section, 
replaces  the  total  field  at  the  surface  by  the  incident  field  times  a  factor  involving 
the  Fresnel  reflection  coefficients.  These  coefficients  are  strictly  valid  only  at 
a  perfectly  smooth,  planar  interface,  but  may  be  employed  with  littlfi  error  when  the 
radii  of  curvature  of  the  surface  at  a  given  point  are  much  larger  than  wavelength. 

When  this  condition  is  met,  the  reflected  field  at  the  given  point  is  treated  as 
though  it  had  been  reflected  from  an  infinite  plane  tangent  to  the  surface  at  that 
point:  hence  the  term  "tangent -plane"  approximation.  The  approximation  is  commonly 
used  to  reduce  an  Integral  equation  for  the  fields  over  the  surface  to  a  definite 
integral,  the  latter  being  called  the  "physical  optics"  integral. 

*The  subject  of  incoherent  scatter  from  a  ground  wave  will  be  discussed  in  a  separate 
report,  and  this  contrast  in  behavior  between  good  and  perfectly  conducting  surfaces 
near  grazing  will  be  examined, 
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la  order  to  be  applicable  for  rough  aurface  acatter  problem* ,  the  radii  of 
curvature  of  the  aurface  at  most  points  must  be  conalderably  greater  than  a  wavelength. 
For  the  ocean,  the  dominant  sea  waves  which  will  contribute  to  HF/VHF  scatter  are 
those  of  the  order  of  one-half  the  radio  wavelength.  To  check  the  applicability,  as¬ 
sume  a  frequency  of  10  MHz,  and  hence,  a  free  space  wavelength  of  30  meters  at  winds 
of  about  10  knots,  15  meter-long  ocean  waves  will  be  excited,  and  their  heights  will 
be  of  the  order  of  C  *  2  to  3  feet.  Such  waves  (assume  they  are  nearly  sinusoidal) 
will  have  a  radius  of  curvature  at  the  peak  and  trough  of  about  7  meters,  considerably 
leas  then  the  radar  wavelength.  Thus,  the  restriction  on  the  radius  of  curvature  at 
HF/VHF  fells,  and  the  use  of  the  tangent-plane  approximation  here  la  extremely 
questionable. 

A  further  argument  against  the  use  of  the  tangent-plane  approximation  with 
the  physical  optics  integral  for  propagation  near  grazing  is  the  inherent  neglect  of 
shadowing  In  the  approximation.  Shadowing  at  these  angles  will  be  a  very  serious 
factor,  and  its  neglect  near  grazing  seems  an  oversimplification. 

Various  investigators  have  nonetheless  employed  this  approximation  and  triad 
to  extend  its  use  of  HF/VHF  scatter  from  the  sea.  As  evidence  of  the  inconsistencies 
resulting  from  its  use,  no  polarization  dependence  is  predicted  for  scatter  within 
the  plane  of  incidence.  Yet  measurements  show  considerably  stronger  scatter  for  the 
vertical  polarization  states  than  for  the  horizontal  states.  For  ground-wave 
propagation,  polarization  is  of  the  essence  of  the  problem,  and  any  theory  which 
shows  no  polarization  dependence  cannot  be  used.  Therefore,  it  is  concluded  thai  the 
tangent-plane  approximation  cannot  be  used  here  for  ground-wave  propagation,  inter¬ 
action,  and  scatter  from  the  sea  in  the  HF/VHF  region. 

4.  Leontovlch  Boundary  Condition 

The  Leontovlch  (or  impedance)  boundary  condition  is  applicable  to  surfaces 
whose  radius  of  curvature  is  larger  than  the  radio  wavelength  inside  the  materials 
as  discussed  in  a  previous  section.  For  the  ocean  this  wavelength  is  »  Xg/I/F^j, 

where  «r  is  the  complex  dielectric  constant  of  sea  water.  At  10  MHz,  this  becomes 
**  1  foot.  Aa  mentioned  in  the  preceding  section,  the  smallest  radius  of  curvature 
for  the  ocean  waves  producing  scatter  at  these  frequencies  is  about  7  meters  or 
21  feet.  Hence,  the  requirement  on  Che  applicability  of  the  Leontovlch  boundary 
condition  is  readily  fulfilled  for  the  HF/VHF  region  and  the  sea  surface. 
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In  employing  this  boundary  condition  along  with  the  Rice  boundary  perturbation 
technique,  one  therefore  avoids  the  necessity  of  writing  the  fields  above  the  surface 
in  terms  of  the  fields  beneath  the  surface.  Since  the  technique  is  now  based  upon 
an  expansion  in  ko£  (l.e.,  the  free-space  radio  wavenumber  times  the  surface  height) 
and  not  in  terms  of  (i.e.,  the  wavenumber  in  sea  water  times  the  surface  height), 
the  difficulties  discussed  in  Subsection  1.  above  are  avoided.  At  HF/VHF,  the 
parameter  koC  is  considerably  less  than  unity  (for  frequencies  below  about  50  MHz  and 
peak-to-trough  wave  heights  below  about  4  feet),  and  the  perturbation  expansion  about 
such  a  parameter  converges  rapidly.  In  addition,  shadowing  is  not  neglected  with  this 
technique,  as  it  is  with  the  tangent  plane  approximation. 

Therefore,  we  conclude  that  the  use  of  the  Leontovich  boundary  condition 
along  with  the  Rice  perturbation  technique  gives  the  soundest  analysis  of  ground-wave 
propagation  and  scatter  across  a  rough  sea  at  HF/VHF. 

C.  Guided  Wave  at  Smooth  Impedance  Boundary 

At  any  smooth,  planar  interface  between  two  media,  various  types  of  aemi- 
infinite  guided  waves  can  be  shown  to  exist.  They  are  semi-infinite  in  the  sense 
that  they  can  extend  to  infinity  in  directions  normal  to  the  interface.  They  represent 
mathematical  solutions  to  Maxwell's  equations  which  are  forced  to  satisfy  the  boundary 
conditions  at  the  interface.  The  existence  of  this  type  of  guided  wave  is  postulated 
without  regard  to  the  reality  of  sources  which  might  be  required  to  excite  them.  In 
this  sense  they  are  analogous  to  the  concept  of  semi-infinite  "plane  waves"  in  free 
space,  the  latter  merely  satisfying  Maxwell's  equations  without  regard  for  the  physical 
reality  of  their  sources.  Like  plane  waves,  the  physical  picture  of  guided  waves  is 
believed  to  be  valid  over  a  small  localized  region  near  an  interface  far  from  the 
sources  producing  them.  As  mentioned  in  Chapter  II,  Wiae^31^  and  Norton^  have 
shown  that  the  Somnerfeld  ground  wave  excited  by  a  finite  source  has  the  same 
polarization  relationship  and  angle  of  incidence  near  the  interface  as  the  guided 
wave. 

The  terra  "guided  wave"  is  employed  here  in  order  to  avoid  the  use  of  the 
controversial  term  "surface  wave".  Actually,  the  guided  waves,  as  we  intend  to  use 
the  term,  is  identical  to  the  "Zenneck  surface  wave"  often  found  in  the  literature. 

Two  types  of  simple  guided  waves  can  be  shown  to  exist  at  a  planar  interface  and 
satisfy  Maxwell's  equations:  plane  guided  waves  and  radial  guided  waves.  The  former 
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represents  the  solution  in  rectangular  (or  Cartesian)  coordinates;  this  wave,  like  a 
plane  wave,  extends  to  Infinity  along  a  line  parallel  to  the  surface  and  perpendicular 
to  the  propagation  direction.  The  radial  wave,  in  contrast,  represents  the  solution 
in  cylindrical  (or  polar)  coordinates;  the  wave  fronts  appear  as  circular,  concentric 
rings  esttnetlng  from  the  z-axis  normal  to  the  Interface.  We  shall  deal  only  with  the 
.former,  or  the  planar  guided  waves  because  of  the  convenience  of  expending  the  surface 
in  rectangular  coordinates. 

Any  elementary  electromagnetics  textbook  can  be  consulted  on  the  derivation 
and  analysis  of  guided  waves.  The  procedure  is  so  standard  that  the  details  will 
not  be  repeated  here.  The  technique  may  be  summarized  aa  follows:  (i)  solutions  to 
the  wave  equation  are  found  in  rectangular  coordinates  for  the  media  on  both  sides  of 
the  interface,  (li)  the  constants  (i.e.,  multiplicative  amplitude  factors  and  also 
wavenumbers  appearing  in  the  exponentials)  are  then  determined  by  applying  the  boundary 
conditions  at  the  interface.  Using  this  technique  and  restricting  attention  to 
Impedance  boundaries  of  Interest  in  this  application,  the  solution  for  the  electric 
field  components  of  the  guided  wave  above  the  surface  can  be  written  as  follows  (see 
Figure  1): 

E^  m  Eo  expliko/l  -  A*  x  -  lk0Az  -  iuit}  ,  (7a) 

E^  “  EodexpUko/l  -  A*  x  -  lk<,Az  -  iu)t)  ,  (7b) 

Hy  I20tt  «*p{iko/l  "  A*  x  -  lk0Az  -  iu»t)  ,  (7c) 


where  Eq  *  E-field  amplitude  constant,  Z#  and  A  are  the  surface 
normalized  surface  impedance  respectively.  The  latter  conforma 
Wait^^’^^,  and  la  defined  as  follows : 


A 


impedance  and 
to  the  notation  of 


(8) 


The  solutions  given  in  (7)  are  based  on  the  assumption  chat  |A|«1  (which  is  true  for 
sea  water  at  HF/VHF).  Note  that  as  the  surface  properties  approach  a  perfect  conductor 
A  0,  and  Equations  (7a,b,c)  degenerate  to  plane-wave  free-spaca  propagation  along  the 
x-dlrection. 


f  12] 

^Jordan's  text  (Section  7.09)  provides  an  analysis  of  guided  waves.  A  more 
detailed  treatment  which  includes  both  planar  and  radial  guided  waves  and  various 
types  of  surfaces  can  be  found  in  Barlow  and  Brown1  (Section  2.1). 
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Figure  1.  Guided  Wave  Above  a  Planar  Impedance  Boundary 

For  a  homogeneous  material  below  the  interface, 


r 


where  and  are  the  relative  permeability  and  permittivity  of  the  medium.  Either 
or  both  may  be  complex.  For  our  application,  sea  water  has  approximately  the  follow¬ 
ing  constants : 

*r  ~  1.  «r  *  80  +  i^  ,  (10) 

where  e0  *  10*^3617,  u>  is  the  angular  frequency,  and  a  is  tha  conductivity  of  aea 
water. *  At  10  MHz,  tfyis  value  is  approximately  *  ■  80  +  17200,  and  thus 
A  *•  1.18  X  10"a  •  e  *4  .  (The- time  dependence,  a  iU>t,  will  be  dropped  henceforth.) 

The  solutions  for  tha  guided  wave  represented  by  Equations  (7a,b,c)  are 
valid  for  a  smooth  sea  surface.  When  the  surface  becomes  rough  or  corrugated,  the 
normalized  surface  impedance.  A,  appears  to  change.  For  a  deterministic  corrugated 
surface,  for  example,  whose  period  is  much  less  than  a  radio  wavelength,  this  change 
manifests  itself  as  a  sharp  increase  in  the  reactive  (or  imaginary)  component  of  A. 

In  particular,  the  surface  appears  inductive.  It  is  this  increase  in  A  which  we  are 

*WaitlJ4J  shows  that  a  more  general  expression  for  the  normalized  surface  impedance 
at  grazing  incidence  above  a  dielectric  surface  is  A  *  l//c  (1  -  1/s  )1/a  ;  this  ex¬ 
pression  is  exact,  and  (9)  is  an  approximation  valid  where  f/|*  |  isrmuch  less  than 
unity,  as  is  true  here.  r 

tThis  conductivity  depends  upon  the  particular  ocean  and  its  temperature.  It  almost 
always  lies  between  3  and  5  mhos/meter.  For  this  report,  a  value  4  mhos/meter  will 
be  used. 
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seeking  for  a  statistically  rough  sea  surface.  This  sverage,  or  effective,  A,  can 
then  be  used  in  the  Norton-Sotnnerfeld  formulation  for  ground-wave  propagation,  in  the 
spirit  of  Wait's  dicusaion^^  . 

D.  Guided  Wave  ac  Slightly  Rough  Impedance  Boundary 

1.  Description  of  Perturbed  Surface 

HflfifeC. 


Let  the  mean  plane  of  the  sea  surface  be  taken  as  the  x-y  plane.  Then  the 

z-coordinate  (or  height)  to  any  surface  point  will  be  designated  £.  See  Figure  2. 

In  general,  £  is  a  function  of  x  and  y.  It  can  now  be  expanded  in  a  two-dimensional 

* 

Fourier  series  over  a  square  area  of  side  length  L.  Thus 

m 

C(x,y)  ■  £  F(«,n)  6xp{ia(mx  +  ny)]  ,  (11) 

s.n*-* 

where  a  ■  —  .  P(m,n)  is  the  Fourier  coefficient  of  the  a,n  th  spatial  harmonic  of 

I*  -«»«=» 

the  surface.  The  real  nature  of  the  surface  height,  C(x,y),  requires  that  P*(m,n)  • 
P(-m,-n),  where  P*  denotes  the  complex  conjugate  of  P.  This  expansion  can  be  employed 
both  for  deterministic  surface-height  profiles,  and  for  random  surface  heights. 


{23] 

*The  notation  of  Rice  will  be  retained  here  as  much  as  possible  in  order  to 
facilitate  reference  to  this  classical  treatment. 
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In  the  case  of  Che  ocean,  £(x,y)  la  a  random  variable.  Over  a  period  of 
Ciae,  Che  surface  profile  undergoes  a  complete  change,  and  we  can  chink  of  an  ensemble 
of  random  height  profiles  valid  aC  different  points  in  time.  Hence,  an  average  over 
time  can  be  physically  justified  as  representing  an  ensemble  average,  without  undue 
discussion  pf  the  ergodic  theorem.  Therefore,  averages  (denoted  here  by  (f)  or  f) 
will  be  understood  to  be  ensemble  averages. 

•  Since  we  shall  have  occasion  to  perform  averages  later,  some  definitions 

will  be  stated.  These  also  facilitate  the  transition  from  a  Fourier  aeries  with 

discrete  surface-height  spectral  components  to  a  continuous,  average  spectrum  of 

[23] 

surface  heights.  The  reader  is  referred  to  Rice1  for  more  detail  on  these 

definitions.  It  can  be  shown  in  any  statistical  treatment  of  random  noise  (e.g.,  aee 
(431  (441 

Davenport  and  Root  or  Rice's  earlier  work  )  that  a  Fourier  series  such  as  (11) 
can  be  used  to  validly  represent  a  random  variable.  Furthermore,  they  show  that  the 
series  coefficients,  P(m,n),  become  uncorrelated  as  L  —  ®.  Practically,  it  is 
necessary  that  L  be  considerably  larger  than  the  surface  correlation  length  or  radius 
in  order  for  the  coefficients  to  appear  uncorrelated;  this  condition  is  assumed 
satisfied  here,  since  the  ocean  area  considered  here  includes  many  ocean  waves. 
Recalling  now  that  the  coordinate  system  was  chosen  so  that  (C(*»y))  “  0,  we  can 
state  that 


<P(m,n)>  -  0 

1*0  for  u,v  ^  -m,-n 


<P(m,n)  P(u,v)> 


W(p,q)  for  u,v  -  -m,-n 


(12a) 


(12b) 


where  p  »  am  ■  2nn/L  and  q  ■  an  ■  2nn/L.  The  function  W(p,q)  defines  the  average 
roughness  spectral  density  of  the  surface,  and  p,q,  are  the  radian  wavenumbers  (or 
spatial  frequencies)  along  the  x-  and  y-dlrections  respectively.  Using  these 
equations,  the  following  relationships  are  established. 


<Ca (*.y)> 


L  <P(m,n)  P(u,v)>  eiax("  +  u>  +  +  v> 

m,n,u,v 


J  r  tt® 

dm  j  dn  r-y  W(p,q) 

m,n  -os  L 
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4  J  I  H<p»q)  dpdq  “  °*  » 


<C(*,y)  C(x\y’>>  -  ;  <?(>».«)  p(u.v)> 

*»n»u,v 


lam  +  iaux*  +  lany  +  lavy' 


-  L  <P(.,n>  P(-»,-n)>  a1*"^  *  X‘>  +  lan(y  ’  y’> 

■7“  .  ..  ? 

cm  m 

"*4J  i  w(p»q>  *ipTx  +  iqTy  dpdq  a  °aR(Tx*Ty>  <l3b 

where  T  ■  %  -  *'  and  T  “  y  -  y 1 .  The  quantity  a®  la  tha  mean-aqua re  height  of  the 

x ,  y 

surface,  and  R(Tx»Tyj  1*  the  surface  height  correlation  coefficient.  Relationship 
(13b)  .aerely  states  that  the  roughness  height  spectral  density  and  surface  height 
correlation  function  are  Fourier  transforms. 


Description  of  Perturbed  Fields 


of  the  Guided  Wave 


The  underlying  philosophy  behind  the  perturbation  approach  to  be  uaed  here 
requires  that  we  expand  the  perturbed  fields  In  the  same  eigenfunctions  as  those  used 
for  the  perturbed  surface.  Such  field  solutions  will  be  required  to  satisfy  tha 
wave  equation.  This  leads  ua  to  chooae  the  following  form  for  our  perturbed  fields: 


E  -  AE(h,0,z)  +  )  A  E(m  +  h,n,z)  , 

x  L-i  mn 

m.n*-1* 


*  -  /  B  E(m  +  h,n,z) 

y  am 

J  m.n*-* 


E  *  E(h,0,z)  +  )  C  E(m  +  h,n,z) 

z  i-t  mn 

where 


E(m  +  h,n,z)  ■  Eo  exp{la(h  +  m)x  +  iany  +  ib(m  +  h,n)z)  , 


BATTSLLS  HIMOMAL  INSTITUTB  -  COLUMBUS  LABORATOMIKS 


25 


b(h  +  m,n)  ■  AJ  -  aa(m  +  h)a  -  a^n*  (13) 

The  definition  of  b  above  is  such  Chat  Equations  (14a, b,c)  satisfy  the  wave 
equation.  The  Cartesian  components  of  Che  H-field  are  not  given  here,  but  are 
readily  determined  from  Maxwell's  equations.  **’" 

In  the  above  equations,  the  presence  of  roughness  manifests  Itself  as  the 

summation  terms.  As  the  roughness  height  approaches  zero,  A  ,  B  ,  C  ,  will  vanish, 

ion  ran  nm 

and  b(h,0)  -  -koA;  Equations  (14a, b,c)  then  become  Identical  with  Equations  (7a,b,c) 
for  a  guided  wave  over  a  smooth  impedance  boundary. 

Also,  Cqq  is  taken  to  be  identically  zero  in  (14c).  This  choice  is  possible; 
what  it  really  means  is  that  all  of  Che  remaining  constants  are  normalized  so  that  the 
0,0  mode  appearing  in  (14c)  is  removed  from  the  summation  as  the  first  term,  E(h,0,z) 
wich  amplitude  Eq.  Physically,  the  guided^ave  portions  of  the  field  appearing 
in  Equations  (14a, b,c)  are  all  terms  having  the  E(h,0,z)  structure.  These 
are 


E®  -  (A  +  AQ0)  E(h,0,z)  ,  (17a) 

-  Bqq  E(h,0,z)  ,  (17b) 

E°  -  E(h,0,z)  .  (17c) 

The  remaining  portions  of  Equations  (14a, b,c)  consist  of  modes  generated 
by  the  roughness.  These  modes,  to  be  termed  the  scattered  field  here,  actually 
include  both  propagating  and  evanescent  modes. 

3,  The  Average  Normalized  Impedance 

The  gulded-vave  portions  of  the  perturbed  fields  are  given  by  Equations 
(17a, b,c).  When  one  compares  (17a)  and  (17c)  with  Equations  (7b)  and  (7a)  for  a  smooth 
Impedance  boundary,  one  is  led  to  define  an  "effective"  or  average  impedance  guided- 
wave  propagation  across  a  roughened  surface  as 

A  »  <A  +  AM)  -  A  +  <A00>  ,  (18) 
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i.e.,  the  affective  impedance  conelste  of  the  constant  Impedance  of  a  smooth  Interface 
plus  (Aqq),  which  accounts  for  the  roughness.  This  convenient  definition  will  in 
fact  be  used,  and  the  goal  o f  the  analysis  will  be  to  derive  an  expression  for  (Aqq). 
If  the  roughness  is  to  contribute  to  the  impedance,  we  expect  this  average  to  be 
nonzero.  Physically,  (18)  says  that  the  average  wave  front  and  polarization  tilt  at 
the  surface  is  (E^)/(£^>  *  A  ■  A  +  (aqq) ,  which  is  another  way  of  defining  the  effec¬ 
tive  surface  impedance. 

Since  there  la  no  E  component  of  a  guided  wave  over  a  smooth  surface,  we 

y  g 

should  expect  that  (B-,.)  appearing  in  (17b)  will  be  zero,  meaning  that  (E  )  is  zero, 

Uv  y 

on  the  average.  This  will  in  fact  be  shown  later  in  the  analysis. 

4,  The  Perturbed  Leohtovlch  Boundary  Condition 


Basically,  the  method  of  solution  to  the  problem  may  be  suouarized  as 

follows:  We  intend  to  employ  Equations  (14a,b,c),  (which  are  solutions  to  the  wave 

equation)  along  with  the  description  of  the  surface  given  in  (11),  to  determine  the 

unknown  constants  A  ,  B  ,  C  ,  and  h,  by  forcing  the  solution  to  satisfy  the 
con  ran  tan 

boundary  condition  at  the  perturbed  surface.  Of  these  constants,  only  is  really 
needed  to  determine  the  effective  surface  Impedance,  as  shown  in  (18).  The  others 
will  prove  to  be  useful,  however,  in  understanding  the  mechanism  of  scatter  from  the 
guided  wave  by  the  roughness. 

The  Leontovlch  boundary  condition,  as  expressed  in  Equation  (3),  is  written 
in  terms  of  the  E  and  H-flelds  at  the  surface,  along  with  the  unit  normal  to  the 
Surface,  n.  All  of  these  quantities  vary  with  position  along  the  surface.  To  reduce 
this  equation  to  a  usable  form,  however,  we  must  express  the  unit  normal  in  terms  of 
Che  surface  slopes;  these  are  the  * rtlal  derivatives  of  the  surface  height  as  a 
function  of  the  Independent  variables,  x  and  y.  Defining  them  as 


r  a  S&teaXl  f  a  ££(*., 

St  dx  ’  *y  dy  * 
this  normal  is  expressed  in  Cartesian  coordinates  as  follows: 


«  «  .  a  a 

-C  x  -  C  y  +  * 

x  *v 

~  *  <i 


(19) 


(20) 
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where  x,  y,  and  z  are  unit  vectors  along  Che  coordinate  axes.  Although  implicit,  and 
are  local  functions  of  position,  x  and  y,  just  as  is  £(x,y). 

With  the  quadratic  function  in  the  denominator,  Equation  (20)  is  not  very 

useful  in  a  perturbation  analysis.  In  order  to  be  useful,  a  further  restriction  on 

the  surface  slopes  must  be  made;  they  must  be  small,  so  that  £*,  Qa  «  1.  When  such 

a  y 

is  the  case,  the  quadratic  can  be  simplified  with  the  binomial  expansion  as  follows 

1 1  +  C*  +  £a]"x/3  “■  l  -  r(Ca  +  Ca),  retaining  only  terms  to  the  second  order. 

Ay  a  x  y 

We  shall  employ  the  above  approximation,  and  substitute  (20)  for  n  into  (3). 

In  addition,  we  shall  express  H  in  terms  of  E  by  using  Maxwell's  equation.  Then 
Equation  (3)  in  vector  form  can  be  reduced  to;  three' scalar  equations  for  the  x,  y,  and 
z  components.  We  preserve  terms  in  £  ,  C  UP  to  the  second  order,  i.e.,  up  to  £*  and 

a  ^  y 

£a.  The  results  for  the  x  and  y  components  are,  respectively., 

.  ,  ,d£  5E  ,  ,•  ,,  ;  .dK  dE 

e.  -  <$*.  -  +  v.  •  f  s,(-if  -  if)  -  L1  -  K? +  «;)]  (if  -  if))  •  <2i*> 

e,  •  v,  -  w. +  v.  -  it  (•  ‘.(if  •  if)  -  -  Ke + o)  iif  -  if))  •  <2ii> 

The  equation  for  the  z  component  is  not  written  above  and  is  not  used  in  the  analysis. 

The  reason  is  that  it  is  not  independent  of  the  other  two,  but  can  be  obtained  from 
them  using  the  divergence  equation. 

Also,  no  attempt  is  made  in  the  above  equations  to  order  the  fields.  For 
instance  E  (the  y-component  at  the  surface)  will  be  small  compared  to  E  ,  i.e.,  it  will 

y  z 

be  at  least  first-order  if  E  is  zero-order.  Hence,  the  term  C  £  E  will  be  third-order 

z  jt*y  y 

and  could  have  been  omitted  from  the  above  equation,  since  only  terms  up  to  and  includ¬ 
ing  second-order  are  to  be  ultimately  retained.  However,  the  field  components  are  not 
ordered  as  to  smallness  at  this  point  so  that  the  possibility  of  omitting  significant 
terma  will  not  occur.  They  can  be  easily  dropped  later,  and  the  intervening  analysis 
will  serve  as  proof  that  they  are  small. 

This  boundary  condition,  as  noted  previously,  applies  at  the  perturbed  surface, 
relating  the  various  E-fleld  components  and  their  derivatives. 
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5.  Restrictions  Made  and  the  Perturbation 
Parameters 


In  many  analyses,  so  many  assumptions  and  approximations  ara  made  along  the 
way  that  the  reader  is  never  sure  of  the  validity  of  the  results.  In  summery,  we 
shall  state  at  this  point  the  assumptions  we  have  made. 


(i)  The  normalized  surface  impedance,  A,  is  considerably  less  than  unity. 
This  is  true  for  sea  water,  and  at  the  highest  frequency  we  Intend  to  consider  (i.e., 
100  MHz),  |Aj  “■  .038;  it  decreases  further  at  lower  frequencies. 

(ii)  The  radii  of  curvature  of  the  ocean  waves  producing  scatter  Is 
considerably  greater  than  the  wavelength  within  the  medium.  At  low  frequenciea,  tills 
restriction  is  the  most  serious,  and  at  10  Milz,  for  example,  It  was  shown  in  a  previous 
section  that  it  is  very  adequately  fulfilled. 

(ill)  In  free  space,  (k0£)a  «  1,  i.e.,  the  surface  height  compared  to 
wavelength  is  small.  In  very  high  seas,  say  Sea  State  6,  where  wave  heights,  £,  can 
be  ±4  feet  from  a  mean  level,  (k0C)3  .4  at  30  MHz.  for  calmer  seas,  which  are  much 

more  prevalent,  the  frequency  of  validity  can  be  extended  higher;  e.g.,  to  100  MHz 
for  See  State  3. 


(iv)  C Ztid  «  I,  i.e.,  the  ocean  waves  responsible  for  scalter  have  rela- 
*  y 

tively  small  slopes.  This  is  adequately  satisfied  for  deep  seas  (i.e.,  away  from 
coastal  reef  areas  which  can  produce  breakers)  where  the  surface  angle  for  waves  whose 
lengths  are  several  meters  will  rarely  exceed  20*,  for  which  <  .14. 


It  is  important  to  note  that  the  above  condltiona  can  even  fai I  at  a  finlLe 
number  of  points  on  the  surface,  and  the  validity  of  the  results  is  not  seriously 
impaired.  They  cannot  fail,  on  the  average,  over  moat  of  the  surface,  however. 

Because  of  the  nature  of  the  perturbation  analysis  which  we  intend  to  employ,  we  will 
have  need  to  "order"  or  expand  terms  mathematically  about  "smallness"  or  "perturbation" 
parameters.  As  a  result  of  the  above  restrictions,  we  shall  employ  the  following 
perturbation  parameters: 


(1)  k0£;  we  shall  include  terms  to  second  order  in  this  |>arameter 

(ii)  Q*  and  ;  likewise,  terms  to  second  order  will  be  retained 

(iii)  A;  we  shall  retain  terms  to  first  order  In  this  parameter 

(lvl  A  .  B  .  C  i  these  are  taken  to  be  of  at  least  first-order  in 
mn  mn  mn  - 
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ama lines*,  a Inc a  they  represent  field  perturbation  due  to  the 

roughness.  On  the  basis  of  the  restrictions  mentioned  previously, 

these  parameters  will  be  shown  to  be  at  least  first-order.  However, 

in  order  to  facilitate  the  derivation,  ^re  shall  employ  the  definition 

A  ■  +  A^  +  ...  ,  with  similar  definitions  for  B  ,  C  .  This 

ntn  mn  on  ion  nm 

means  that  we  break  A  up  into  firat-order  contributions,  second-order 
mn 

contributions,  and  so  on.  The  first-order  contributions  to  A 

wa 

(i.e.,  A^  )  then  come  from  first-order  terms  in  k0C,  Cy»  and  A. 

Derivation  of  First-Order  Coefficients 

of  the  Perturbed  Field:  a(*  \  B^),  and  c(i) 

_  am  mn  _ am 


Equations  (14a,b,c)  represent  the  perturbed  field  expressed  in  terms  of  the 

coefficients  A  ,  B  ,  and  C  .  These  coefficients  we  order  according  to  smallness, 
mn  mn  on 

as  discussed  in  the  preceding  paragraph.  Then  the  resulting  expressions  for  the 
perturbed  fields  are  substituted  into  the  boundary  conditions  (21)  and  evaluated  at 
the  boundary,  i.e.,  for  z  ■  £.  At  this  point,  we  also  discard  terms  obviously  of 
higher  order  than  second.  We  employ  the  following  expansions  of  the  fields  at  the 
surface,  correct  to  the  second  order  in  £: 


E(m  +  h,n,z) 


E(m  +  h,n,C)  *  F(h)  Ex(m,n)  [l  +  ib(m  +  h,n)C  -  |  b»  (m  +  h,n)Ca]  ,  (22a) 


3E(m  +  h.n.z) 

3x  z 


i(h  +  m)a  F(h)  Ex(m,n)  £l  +  ib(m  +  h,n)C  -  j  i?  (m  +  h,n)Ca  , 


dE(m  +  h.n.z) 

9y  z-C 


ina  F(h)  Ex(m,n)  £l  +  ib(m+h,n)C  -  j  b*  (m  +  h,n)£a  j  , 
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lb(m  +  h,n)  F(h)  Ex(ni,n)  [l  +  lb(ta  +  h,n){  -  j  &  (m  +  h,n)C*]  ,  (22d) 

where  F(h)  *  EqC  ,  end  Ex(n,n)  "  axpfiamx  +  iany ] . 

When  Che  above  expansions  are  used,  the  follcving  expressions  are  obtained 
from  (21),  correct  to  the  second  order  In  the  smallness  parameters  defined  In  the 
preceding  section. 


A  F(h)  [l  +  lb(h,0)C  -  |  b*(h,0)Ca]  + 

+  X{[A«  +  A«  ][l  ■  “><"  +  h,n)C  -  7^(m  +  h,n)C*]}  . 

•  F(h)  Ex(m,n)  -  Q  ?0»)  +  C/0»>  [l  +  lb(h,0)c]  +  C,  ^  Cm>  F(h)  Ex<a»n>  " 

"  [<y  i  iB»>  1(m  +  h)*  ’  in-}  F(h)  Ex<n*n)  ’ 

-  ib(h.O)  A  F(h)  [l  +  lb(h,0)C  -  f  b»(h,0)C*]  + 

+  lha  F(h)  [l  +  lb(h.0)C  -  |  ba(h,0)C*  j  -  £  {  [  (A^)  +  A^})  ib(m  +  h,n)  - 


*  (C«  +  Cm)  1("  +  h>«]  ’ 

•  [l  +  ib(m  +  h,n)c]  }  F(h)  Ex(m,n)  +  j  +  Cj)  [ib(h,0)A  -  lha]  F(h)]  ,  (23a) 
1  {  [Bli>  +  [l  +  !*»(»  +  h,n)c]  }  F(h)  E*(m,n)  -  CXC/  F(h)  + 

+  Cy  F(h)[l  +  ib(h,0)cj 

+  Cy  XC™  F(h)  Ex<-»“>  “  it  L"Cx  I  {  lM  “  Bli>  *<■  +  h>‘)  F<h>  *»<■.“> 

-  L  {  [*£  +  B™)  ib<B  +  h*n)  *  (C™  +  C^)  H  [l  +  ib<“  +  h'")C]  } 

F(h)  Ex(tn,n)J  .  (23b) 
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la  tha  above  equations,  the  factor  F(h)  ia  common  to  all  terms  and  can  be  dropped. 

Following  the  classical  perturbation  technique,  we  now  collect  terms  of  each 
order  and  set  each  set  separately  equal  to  zero.  Let  us  collect  first  the  zero-order 
terms  from  (23a)  (there  are  no  zero-order  terms  in  (23b)). 

A  +  ~  [ib(h,0)A  -  iha]  -  0  (24) 

This  aquation  contains  the  propagation  constants  associated  with  the  guided  wave.  In 
order  to  convince  ourselves  that  the  left  side  is  really  equal  to  zero  (to  at  least 
the  second  order),  let  us  employ  the  divergence  condition  on  the  gulded-wave  fields 
as  expressed  in  (17). 

—  3E  G  oE  G  3E  G 

7.eG.0-_2L_  +  _J£_  +  _S_.  lha(A  +  Aqo)  +  ib(h,0)  .  o  , 


b(h,0)  -  -ha (A  +  A^j) 


(25) 


But,  ha  is  defined  in  terms  of  b(h,0)  from  (16),  i.e.,  ha  -  /ko  -  ba(h,0).  The 
quantity  b(h,0)  will  be  small  compared  to  ko  since  it  represents  the  z-portion  of 

* 

the  wavenumber  for  the  guided  wave.  On  the  other  hand,  ha  will  be  very  close  to  ko  . 
Hence 


ha 


Ml 


1  b?jh.QJ, 

2  koa  +  ***' 


(26) 


by  the  binomial  expansion.  How,  let  us  substitute  (25)  for  b(h,0)  into  (24). 

a  -  i£ 1 IV  ♦  W + 1]  -»  • 


New,  substitute  (26)  for  ha  into  the  above.  The  result  is 


A* 


+  A*  A 


00 


1  ba(h .0) 

2  ko3 


A 


0 


(27) 


All  of  the  terms  in  the  above  equation  are  of  higher  order  than  second  in  A,  A^,  and 
b(h,0)  .  Hence,  Equation  (24)  for  the  zero-order  terms  is  actually  zero  to  the  order  of 
terms  to  be  retained  throughout  this  analysis. 


*These  statements  are  true  for  the  types  of  surfaces  being  considered  here,  i.e.,  those 
which  satisfy  the  restrictions  imposed  in  the  preceding  section. 
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Mow,  lac  u*  collect  all  of  cha  tarma  from  (23a)  and  (23b)  which  ara  of 
firat-order  in  k<,C.  Cy>  «nd  C^.  Thay  ara: 

1A  b(h,0)C  +  C*  +  £  iba(h,0)C  -  iha  b(h,0)C  + 1  {  [l  +  b<m  *  b.«»>]  * 

-  A  (m  +  h)a  fW)  Ex(m,n)  -  0  ,  (28a) 


<y + 1  { b +  £ b(m  +  h*n>]  B«>  *  £  “  C«1  Ex<B'n> " 0  •  (28b) 

In  aquation  (28a)  above,  Cho  terms  multiplying  ib(h,0)£  are  identically 
the  same  terma  aa  appeared  in  Equation  (24);  they  are  equal  to  zero,  at  least  out 
beyond  the  second  order.  Henca,  Equation  (28a)  becomes 

t.tI{[1  +  ^,,<“th-">]A™  (■  +  ».  c™}  »(.,»>- o  •  <“•=> 


The  surface  slopes  C  and  £  can  be  written  as  a  series  in  tha  eigenfunctions 
tc  y 

Ex(m,n)  by  differentiating  Equation  (11): 


ima  P(m,n)  Ex(m,n); 


l 


ina  P(m,n)  Ex(m,n) 


(29) 


When  these  expressions  are  aubatituted  into  (28b)  and  (28c),  one  has  a 

sariaa  whose  terms  each  contain  tha  same  eigenfunction  Ex(m,n)  aa  a  factor.  In  order 

* 

for  such  a  sariaa  to  be  zero,  each  term  must  be  identically  equal  to  zero  .  Tha 
following  two  equations  are  obtained  from  this  procedure. 


+  &L±&£l]  a<1>  .  a(B_±, cO>  .  .lM  H  a) 

L  Kq  -l  mn  Kq  am 


(30e) 


fl  +  Ab<-W  -  B(l) 

L  ko  mn 


,U> 


-ins  P(m,n) 


(30b) 


Represented  above  are  two  equations  in  three  unknowns,  and 

rev  no  am 

C'  .  To  this  we  can  add  a  third  Independent  equation  by  using  the  divergence 
con  _ 

condition  (i.a.,  V  •  E  ■  0)  with  Equations  (14e,b,c).  Tha  m,n  *  0  terma  give  (25). 
All  higher  terms  when  equated  to  zero  give  the  following  identity: 


,(l> 


(m  +  h)  a  Av  '  +  na  Bv  '  +  b(m  +  h.n)  C 


(1) 


mn 


,(1) 

mn 


(30c) 


*This  result  can  be  proved  formally  by  integrating  the  series  over  the  range  of  x  and 
y,  i.a.,  from  -L/2  *  x,y  *  L/2. 
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Now  we  have  three  equations  in  three  unknowns.  They  are  most  easily  solved 

by  expressing  in  terms  of  A^  and  B^  using  (30c).  This  then  when  substituted 

ran  ran  ran 

into  (30a)  and  (30b)  leaves  two  equations  in  two  unknowns,  which  are  then  readily 
solved  by  elimination.  The  final  results  for  the  first-order  perturbed  field  coef- 
ficients  are: 


.<!>  .  Jh. 


D  P(m’n);  Bmn 


(1)  -Tp(m«n>;  C« 


-(m+h)aN.  -  naN 
A  D 

b(m  +  h,n)D 


P(m,n) 


where 


-  -H) -  • 

» •  0  ♦  ‘PHM1  +  ££  ?&)]  [> ♦  ‘PHM1  *  • 


(31) 

(32a) 

(32b) 


A^(m  +  h)an3a4 
k£bJ(m  +  h,n) 


(32c) 


Several  facts  are  worth  noting  in  Che  above  equations  for  the  first-order 
coefficients  of  the  perturbed  fields.  First  of  all,  the  quantities  N  ,  N  ,  and  D 
contained  in  the  factor  before  P(m,n)  are  functions  of  the  scattered-field  propagation 

directions  and  Che  surface  material  properties;  the  former  directions  are  determined 

J 

by  (m  +  h)a  and  na,  while  Che  surface  me  ter la 1  properties  are  represented  entirely 

by  the  normalized  surface  impedance,  A  “  J  .  These  quantities  do  not  depend  upon 

the  surface  roughness  profile.  The  dependence  upon  this  profile  is  contained  entirely 

in  the  factor  P(m,n).  The  latter,  recall,  is  the  m,n-th  Fourier  coefficient  in  the 

expansion  for  the  surface  height,  as  expressed  in  (11).  Hence,  the  strength  of  the 

scattered  plane  wave  in  Che  direction  defined  by  (a  +  h)a,  na,  is  directly  proportional 

to  the  m,n-th  Fourier  component  of  the  surface  height  (to  the  first  order).  If  a 

surface-height  component  corresponding  to  m,n  is  zero,  then  A^l\  B^,  and  are 

ran  ran  ran 

zero,  meaning  Chat  no  scattered  wave  propagates  in  the  direction  defined  by  (m  +  h)a,na. 

The  theory  of  scatter  from  slightly-rough  periodic  surface  dates  back  to  Lord 
[45] 

Rayleigh  .  He  derived  expressions  similar  to  (31)  and  noted  that  scatter  took 
place  in  unique  directions,  or  lobes,  determined  by  the  Fourier  components  of  the 
surface  and  the  length  of  their  periods  with  respect  to  the  wavelength. 
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Lat  us  recall  chat  Che  object  o f  Che  present  analysis  is  the  derivation  of 

aqqJ  ^  the  average  of  this  coefficient  which  contributes  to  the  surface  Impedance, 

as  seen  in  (18).  To  the  first  order,  however,  is  zero.  This  is  true  because 

A00^  ia  directly  proportional  to  P(0,0),  i.a.,  the  DC  or  naan  height  of  the  surface 

C(*»y>  above  the  x-y  plane.  The  x-y  plane  was  chosen  to  be  the  naan  plane,  however, 

and  hence  P(0,0)  is  identically  zero.  Thus  we  see  that  A  can  be  no  greater  than 

00 

second  order.  It  is  for  this  reason  that  second-order  terms  were  retained  in  the 
equations. 

7.  Derivation  of  the  Second-Order  Coefficient  A*2* 


As  shown  in  the  preceding  section,  the  contribution  of  the  roughness  to 
the  surface  impedance  is  contained  in  A*2);  A ^  is  identically  zero.  Hence,  we 
shall  determine  A^  in  this  section  from  Equations  (23).  While  it  is  possible  to 
determine  all  of  the  general  second-order  coefficients,  A^2\  c^2\  this  will 

not  be  done  here;  they  are  not  necessary  for  our  purposes. 

(2) 

It  turns  out  that  AM'  can  be  determined  entirely  from  (23a).  Keeping  with 
the  perturbation  technique,  let  us  gather  the  second-order  terms  in  (23a)  and  equate 
them  to  zero. 

-  J  »a(h,0)Ca[4  -  £  b(h.O)  +  £  ha]  +  l  [ib(«  +  h,n)CA^>  +  A^>]  Ex(.,n) 

‘  <  +  lb<h-°>«*  +  Ic™  «*<■.■>  "  £  [-  C yL  W”  iC  ♦  h>. 

A™>  iM}  Ex<m’n)  -  I  {lb<”  +  h*n)c[A^)ib(m  +  h,n)  -  <£>  i(m  +  h)n]  +  [a^IMm  +  h,n) 

"  +  h>n]  }  Ex(«,n)  +  £  (fij  +  Cj)  (ib(h,0)4  -  iha)J  .  (33] 

The  first  term  containing  £*  is  zero  because  the  factor  in  square  brackets 
is  zero,  as  seen  from  (24).  Also,  using  (24),  the  last  term  can  be  simplified,  since 
the  factor  ib(h,0)A  -  iha  is  -ik<,. 

Now,  let  us  group  the  unknown,  second-order  coefficients  A^2*  and  C*2* 

mil  mu 

appearing  in  (33)  together  on  the  left  side  of  the  equation.  The  result  is 
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IK*-*  6iat!usl)  • »*£*<£>]  “<».»>  •}(«-*)♦  wn 

- 1  {  [ib(.  +  h,„>c  (i  ♦  4  *J!Lt±Al) -  4  tt  g  Aa>  „  4  laiSJi  Cy  ,W 

+  [-lb<»  +  h,n)  CA  ^  +  tj  <£>}  E*(m,n)  .  ( 

The  right  aide  of  the  above  equation  will  involve  double  aummation  seta. 
In  order  to  employ  the  orthogonality  relationship  of  the  eigenfunctions,  Ex(m,n), 
more  effectively  it  ia  convenient  to  rearrange  the  double  summation.  The  examples 
below  illustrate  the  procedure. 

C!  -  7  7  (i«0  (iVh)  P(cr,P)  P(Y,C)  Ex(o  +  Y,  5  +  0 


^  (i<»)  (iVa)  P(cr,0)  P(Y,0  &(<*  +  Y,  0  +  0 

a  y  )  la(m  -  A)  (iaA)  P(A,1)  P(m  -  A,  n  -  i)  Ex(m,n)  ,  (35a) 

m,n  M 

7  ib(m  +  h.nJCA^  -77  ib(a  +  h,  0)  P(Y,0  A<J>  Ex(or  +  Y.  0  +  0 
m,n  YTfc 

*  y  y  ib(A  +  h,l)  P(m  -  A,  n  -  1)  Ex(m,n) 
m,n 


ib(A  +  h,l)  P(m  -  A,  n  -  1)  Ex(m,n) 


Likewise,  the  remainder  of  the  terms  on  the  right  side  can  be  arranged  so  that  the 
factor  Ex(m,n)  appears  expliclty.  Aa  a  result,  the  right  side  of  (34)  then  becomea: 

|77[(iaA)  ia(m  -  A)  -  (iaf)  ia(n  -  i)]  P(A,Z)  p(n  -  k.n  -  i)  Ex(m,n) 

+  ikoA  77  i*(m  “  *)  P(A, t)  P(m  -  A,  n  -  i)  Ex(m,n) 

-  77{a^  [ib(A  +  h,J)  +  A  A  ■■  <1  ]  +  B(l)  A  iaaCA^h}(n.-.fl 

+  cj^  [ia(m  -  k)  -  A  lbl*±  MM*  +  h)]}  P(ra  -  A,n  -  1)  Ex(m,n)  .  (36) 
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Recall,  we  ere  interested  in  finding  A ^  ,  i.e.,  the  second-order  coefficient 
for  a  ■  n  ■  0.  This  can  be  readily  done  by  invoking  the  orthogonality  of  the  Ex(m,n) 
functions  over  the  square  region  of  interest.  As  a  result  of  this  relationship,  the 
coefficients  of  Ex(0,0)  on  each  side  of  che  aquation  are  equal. 


The  left  side  of  (34)  becomes 


»«>  .  4  6i  c<?> 

00  ko  00 


This  result  can  be  simplified  further.  First  of  all,  Cqq  is  zero,  by  definition,  to 

all  orders;  this  definition  is  merely  one  possible  normalization,  as  discussed  after 

Equations  (14).  As  seen  above,  it  Is  a  logical  selection  because  it  results  in  a 

convenient  separation  of  coefficients.  Secondly,  the  term  A  is  of  second  order 

compered  to  unity.  This  is  evident  from  (23)  where  we  sea  that  =- 

“  *0 

“>-An ■ 


Hence,  it  can  be  neglected  compared  to  unity,  to  the  order  of  the  analysis 


here.  Then  the  left  side  becomes  merely  A 


(2) 

00 


When  we  group  the  terms  on  the  right  side  of  (34),  as  shown  in  (36),  which 


multiply  Ex(0,0),  we  have  the  resulting  equation  for  A 


(2) 

00 


*0?  -  I  .  {[-ib<*  +  M>  +  cu]  +  *  ♦  —  i— ]  H\l) 


A  [(*££  +  A<i>  +  iaZULtm  ,<J)  +  ML +  c<i >]}*<-*> 

In  the  above  equation,  the  terms  in  the  first  set  of  square^ brackets 


(37) 


within  the  summation  are  the  lowest-order  terms  in  the  surface^  Impedance,  A.  For  a 
perfectly  conducting  ro'igh  surface,  A  0,  and  only  these  first  terms  are  left.  Hence, 
the  remaining  terms  represent  the  contribution  due  to  the  finite  conductivity  of  the 
surface  material,  correct  to  che  first  order  in  A. 

At  this  point,  we  can  employ  the  expressions  derived  earlier  for  the  first- 

order  coefficients  A^\  and  C^;  these  are  shown  In  Equations  (31)  and  (32). 

mn  mn  ton 

When  they  are  substituted  into  (37),  considerable  algebraic  simplification  is  possible. 
Sperlng  the  reader  the  details,  we  write  the  final  expressions  for  tha  various  terms. 

In  reducing  tha  expression  in  the  first  set  of  square  brackets,  we  preserve 
terms  up  to  A.  The  final  result  is 
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tb(A  +  h,l)  +  iaA  C 


,(D 

'Af 


L- 


kpaa  A9 


(A  +  h,i)D(A, i)  T  D(A, t) 


(aaA*  +  a9fa)]  P( A, i) 


(38a) 


where  D(A,Z)  is  given  in  (32c)  with  m,n  there  replaced  by  Ati. 

The  expression  in  the  second  set  of  square  brackets  is  already  expressed 
in  the  form  we  desire.  That  in  the  third  set  of  brackets  can  be  simplified,  and  only 
terms  of  zero-order  in  A  need  be  retained  because  the  expression  la  already  of  order 
A  dud  to  this  multiplicative  factor.  Thus  to  the  lowest  order,  the  expression  within 
these  third  brackets  becomes 

[5^17+  0(A)]  P(M)  .  (38b) 

Using  the  simplifications  expressed  in  (38),  Equation  (37)  can  thus  be 
written  as  follows: 

.(2)  »  V  J  kQaaAa  fa^  +  aaia  -  k»aA 

*00  fa  \b(A  +  h,f)D( A, /)  +  a  L  D(A, i) 


+  (koaA  +  aa-ia)jj  |P( A, f)  |a  .  (39) 

In  addition,  the  expression  D(A,f)  can  be  reduced,  1  nd  to  order  A,  it 

becomes 


D(A,f) 


1  + 


1*0 

b( A  +  h,f) 


(AO) 


(2) 

Thus,  we  have  obtained  an  expression  for  Aqq  in  terms  of  a  summation  of 
terma  in  the  square  of  the  Fourier  coefficients  of  the  surface  height.  This  constant, 
as  shown  in  (18),  is  the  contribution  to  the  surface  impedance  due  to  the  presence 
of  the  roughness.  Since  terma  are  retained  to  order  A  in  aqq\  one  can  observe  the 
effects  of  finite  conductivity  of  the  surface  material  on  the  roughness  contribution. 
When  the  material  is  such  that  |a|  Is  large,  Equation  (39)  is  obviously  Inadequate 
to  describe  the  effect  of  roughness.  The  second  term  in  (39)  provides  a  first-order 
correction  in  A  for  moderately  small  values  of  this  parameter. 
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Equation  (18)  expresses  the  overall  effective  surface  impedance  aa  the 

sum  of  A,  the  constant  impedance  of  the  surface  material  in  the  absence  of  roughness, 
(21 

pitta  the  average  of  aqq  »  which  is  the  effective  contribution  due  to  roughness. 

Before  proceeding  to  Cake  the  indicated  average,  it  should  be  noted  Chat  there  may 

(21 

be  some  cases  where  an  average  is  not  desired.  Equation  (39)  expresses  aqq  exactly, 

in  terms  the  coefficients  of  the  Fourier  expansion  for  the  surface  profile.  If  this 

surface  height  profile  is  known  exactly,  these  coefficients  can  be  determined  exactly, 

* 

and  (39)  represents  a  deterministic  result  •  For  example,  if  the  surface  profile  is 

a  pure,  one-dimensional  sinusoid,  then  £(x)  can  be  expanded  so  that  only  two  terms  In 

(2) 

the  series  (11)  are  non-zero.  In  this  case,  consists  of  only  two  terms,  the  one 
which  has  ft  •  ±1,  i  ■  0. 

/ 

The  sea  surface  is  an  example  in  which  Che  exact  height  profile  is  neither 
known  deterministically  nor  constant  over  a  very  long  period  of  time.  Hence,  it  is 
one  situation  in  which  an  average  is  the  only  meaningful  description  of  the  effects 
of  the  surface  on  propagation.  To  take  the  average  of  (39),  we  shall  employ  the 
definition  of  the  height  spectrun  W(p,q)  in  terms  of  (|p(m,n)|a),  as  expressed  In 
(12b).  In  addition,  we  allow  the  animation  over  A,1  to  approach  an  Integral,  in  the 
same  manner  aa  was  done  in  (13).  Understanding  here  chat  p  ■  aft  -  2nA/L,  and 
q  -  al  «  2rri/L,  we  obtain  the  following  result  for  <A^q>,  and  hence  A: 


A 


A  +  A, 


(2) 

00 


* 


or 


A  *  A  + 


Iff 


*(P.9)  W(p,q)dpdq 


where 


?(„  q)  .  E?-±  ZIMeL+JL >o£)  +  _ 

b'  +  A(b'a  +  1)  2 


(41) 


(42a) 


b'  -  r-Ao  -  (p  +  k0)"  -  qJ  ,  (42b) 


■*The  validity  of  the  result  is  still  dependent  upon  fulfillment  of  the  restrictions 
described  In  Section  5  above,  of  course. 
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and  W(p,q)  ia  Che  two-dimens tonal  surface  height  roughness  spectrum;  it  is  a  function 
the  spatial  wavenumbers  p,q  corresponding  to  Che  x,y  directions  and  is  defined  in 
terms  of  the  height  correlation  function  in  (13b),  In  obtaining  the  above  expressions, 
we  also  used  the  approximation  that  ha  -*  kg;  this  is  valid  so  long  as  4  is  small 
compared  to  unity,  as  seen  from  (26). 

As  a  check  on  (41),  if  we  permit  4  to  approach  zero  in  (42a),  we  obtain 
F(p,q)  ■  Pa/b',  and 


i  rVfi 

4  J  J 


this  result  was  derived  from  an  integral  equation  technique  by  Feinberg  in  19441  , 

and  repeated  by  Brenner1  1 .  In  addition,  if  we  permit  the  surface  to  become  perfectly 

conducting  so  that  4  approaches  zero  everywhere,  (43)  consists  only  of  the  integral; 

r  23] 

this  expression  checks  with  that  of  Rice  ,  who  Initially  assumed  a  perfectly 
conducting  surface.  Hence,  we  have  in  (41)  an  expression  valid  to  order  4  for 
imperfectly  conducting,  rough  surfaces. 

9.  Physical  Interpretation  of  Roughness  Contribution 
to  Surface  Impedance 


Equations  (41)  and  (43)  show  that  the  contribution  of  the  roughness  to  the 
surface  impedance  la  represented  by  the  integral.  To  better  understand  the  inter¬ 
action  of  the  guided  wave  with  the  roughness,  let  us  consider  first  the  simpler  form 
of  the  integral,  l.e.,  that  in  (43),  which  is  valid  when  the  surface  is  perfectly 
conducting.  The  height  spectrisn,  W(p,q),  ia  always  a  positive  real  quantity.  Thus 
the  nature  of  the  Integra ^contribution  depends  entirely  upon  the  denominator,  i.a., 
b'  ■  Jl  -  +  l^3  -  .  This  quantity  can  be  either  gure  real  and  positive  or 

purely  imaginary  and  positive,  depending  upon  whether  ia  l*88  than 

or  greater  than  unity. 


This  can  be  better  illustrated  by  referring  to  Figure  3,  Let  us  separate 
the  contribution  to  4  of  the  Integral  into  two  parts,  R^  -  1X^.  The  contribution  R^ 
comes  entirely  from  that  part  of  the  height  spectrum  lying  within  the  unit  circle 
centered  it  ^-  *  -1,  “•  ■  0.  If  there  are  no  roughness  spatial  frequencies,  or  waves, 

Kq  Kg 

within  this  region,  then  the  resistive  contribution  is  zero;  these  spatial  frequencies 


■ATTELLE  MEMORIAL  INSTITUTE  -  COLUMBUS  LABORATORIES 


Spectral  region  in  which  surface  roughness  produces 


BATTELLE  MEMORIAL  INSTI TUTE  -  COLUMBUS  LABORATORIES 


best  AVAILABLE  copv 


Figure  3.  Effect  of  Various  Regions  of  Spatial  Roughness 
Spectrum  on  Effective  Surface  Impedance 
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must  thus  b«  less  than  2ko.  Physically,  this  statement  means  that  only  ocean  waves 

longer  than  r  can  contribute  to  the  resistive  portion.  From  a  study  of  scatter  from 
1  [41  42] 

this  type  of  surface  ’  ,  it  has  been  shown  that  waves  whose  lengths  are  greater 

than  X/2  are  responsible  for  scatter.  Hence  the  interpretation  of  the  resistive 
portion  becomes  clearer;  longer  ocean  waves  whose  wavenumbers  lie  within  the  circle 
are  responsible  for  removal  of  energy  from  the  guided  wave  and  scatter  of  this  energy 
into  all  directions  in  the  upper  hemisphere.  This  energy  removal  produces  an  increase 
in  the  resistive  term  of  the  surface  impedance. 

On  the  other  hand,  if  there  are  no  ocean  waves  whose  lengths  are  greater 

than  a  ha  If '-wave  length  (i.e.,  which  lie  outside  the  unit  circle),  then  the  roughness 

contribution  to  A  is  purely  reactive;  in  addition,  it  is  always  an  inductive 

[41  421 

reactance.  From  scatter  theory  *  ,  roughness  waves  of  these  shorter  lengths  do 

not  scatter  (at  least  to  the  first  order  in  k<,C,  the  roughness  height).  Hence,  this 

higher  frequency  roughness  produces  a  perturbation  on  the  local  field  at  the  surface 

which  exists  only  at  and  near  the  region  between  the  waves;  since  there  are  no 

scattered  propagating  fields  removing  energy  from  the  guided  wave,  this  effect  should 

be  evident  only  very  near  the  surface.  The  perturbed  modes  in  this  case  are  not 

propagating,  but  evanescent,  because  the  only  coefficients  and 

mn  mn  ’  mn 

which  are  non-zero  are  those  whose  wavenumbers  am,  an  are  such  that  b(m  +  h,n)  is 
imaginary,  and  hence,  the  exponential  e*'^®  +  *1,n^z  in  the  perturbed  fields  (Equations 
14-16)  is  attenuating  in  the  +z  direction  above  the  surface. 

This  reactive  nature  of  surfaces  with  short  roughness  periods  is  confirmed 

by  many  studies  on  corrugated  surfaces,  and  especially,  surfaces  with  rectangular 

slots  or  rlbs^  ’  Elliott's  analysia^^  is  baaed  upon  Floquet's  theorem  (as  are 

many  other  similar  treatments),  and  employs  an  approximation  of  the  field  between 

the  ribs.  The  approximation  is  valid,  Elliott  points  out,  when  there  are  more  than 

f  331 

ten  corrugation  periods  per  radio  wavelength.  Wait1  1  feels  this  requirement  can 
be  relaxed  to  five  periods  per  wavelength.  These  analyses  then  show  that  the  dominant 
guided  mode  has  a  wavenumber  which  can  be  related  to  an  effective,  purely  reactive 
(inductive),  surface  impedance.  The  slots  or  ribs  act  as  resonators  in  this  case. 

The  analysis  we  used  here  cannot  be  used  to  analyze  such  structures,  because  of  the 
requirements  we  imposed  in  Section  5  above  demards  that  the  surface  slopes  be  small. 

A  ribbed  surface  has  Infinite  surface  slopes  at  the  edges  or  walls  of  the  ribs. 
Nonetheless,  the  two  different  techniques  both  lead  to  reactive  surface  impedances 
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when  there  are  more  then  two  surface  periods  per  each  radio  wavelength.  We  shell  see 
later  that  the  reactances  predicted  by  our  analysis  are  well  belcw  those  applicable 
to  ribbed  surfaces  with  siaiilar  heights  but  whose  walls  sre  vertical;  such  a  result 
is  expected,  because  more  sloping  sides  will  reduce  the  resonator  offset  between 
teeth. 

When  the  surface  is  not  perfectly  conducting,  finite  fields  can  exist  beneath 
the  surface  as  well  as  above  it.  This  fact  is  accounted  for  by  the  Leontovich 
boundary  condition  and  the  use  of  a  non-zero  surface  impedance,  A^A  ■  J  \  One  ef¬ 
fect  of  this  finite  impedancs  is  in  evidence  in  the  denominator  of  the  Integra  in  (41) 
the  denominator  does  not  go  to  zero  on  the  unit  circle  in  Figure  3.  This  fact  nukes 
for  easier  numerical  integration  in  the  vicinity  of  this  singularity.  Results  later 
will  show,  however,  that  the  effect  of  the  additional  terms  in  (42a)  is  to  increase 
slightly  both  the  resistive  and  reactive  portions  of  the  effective  surface  impedance. 


Numerical  Determination  of  Effective  Surface  Impedance 
for  Two  Ocean  Wind-Wave  Models 


The  roughness  present  on  the  ocean  surface  is  due  almost  entirely  to  the 
wind  blowing  across  the  water.  Oceanographers  have  found  that  to  excite  ocean  waves 
of  a  given  spatial  period,  L,  a  wind  with  velocity  greater  than  J  ^  muat  blow  for 
several  hours  (g  is  the  acceleration  of  gravity).  If  the  wind  does  not  reach  this 
given  speed,  then  waves  of  this  length  are  simply  not  present.  If  the  wind  speed  doea 
exceed  this  speed,  the  heights  of  the  waves  of  period  L  remain  constant,  relatively 
unaffected  by  any  further  wind  increase.  Hence  the  shorter  ocean  waves  are  the  first 
to  be  excited  by  an  increasing  wind,  but  after  aroused,  the  height  of  these  short 
waves  remains  fixed.  In  this  sense  the  ocean  wave-height  spectrum  is  said  to  saturate 
with  wind  speed. 

The  saturated  height  vs  spatial  period  for  ocean  waves  has  been  the  subject 

of  experimental  and  theoretical  studies  to  oceanographers  for  many  years.  Two  recent 

monographs  on  ocean  waves  provide  an  excellent  discussion  and  review  of  experimental 

35  47 

evidence  on  thia  subject  (Kinsmen  and  Phillips  ).  Over  the  pest  15  years,  several 
empirical  models  have  been  proposed,  and  two  of  these  have  found  relatively  wide 
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acceptance'.  They  differ  mainly  in  the  "power  law"  followed  at  aaturation  and  in  the 
ahape  of  the  longwave  cutoff  vs  wind  speed.  They  are  the  Neumann-Pierson  and  the 
Phillips  spectra. 

It  is  not  entirely  correct  to  say  that  the  ocean-wave  heights  at  a  given 
point  on  the  sea  are  due  solely  to  the  winds  that  have  blown  across  this  point. 

Waves  generated  in  a  given  area  can  and  do  propagate  to  portions  of  the  ocean  which 
have  had  considerably  different  wind  histories.  This  is  especially  true  of  the  longer, 
faster  moving  ocean  waves.  The  roughness  in  a  given  area  due  to  winds  in  other  areas 
is  known  as  swell.  The  ocean-wave  spectra  models  to  be  discussed  below  neglect  swell 
as  a  source  of  roughness,  and  relate  the  ocean-wave  heights  to  the  winda  above  that 
surface  area.  This  oversimplification  should  be  kept  in  mind  before  placing  too  much 
confidence  in  the  wind-dependence  of  the  effective  surface  Impedance. 

2.  The  Neumann-Pierson  Spectrum 

Neumann  and  Pierson  first  suggested  a  model  which  seemed  to  fit  wavemeter 

★ 

measurements  about  IS  years  ago.  Their  model  was  an  isotropic  temporal  spectrum  in 

which  the  saturated  value  is  related  to  the  temporal  radian  wave  frequency  as  u>~*. 

To  this  they  added  a  "guesstimated"  lower  end  continuous  cutoff  function,  since  it 

seemed  unreasonable  to  them  that  the  cutoff  should  be  abrupt.  The  resulting  empirical 

*  35 

law  can  be  converted  to  an  isotropic  spatial  spectrum.  Kinsman  proposes  a  reason- 

if 

able  extension  of  this  to  produce  a  directional  spatial  spectrum  with  a  cosine- 
squared  dependence  about  the  wind  direction  in  azimuth.  The  spatial  spectrum  then 
has  the  form 


w(ptq) 


C(p  coscr  +  q  sinor)a 


«p{-2g/(UVpSr+'qT)  j 


(44) 


where  U  is  the  wind  speed,  g  is  the  acceleration  of  gravity  (9.81  m/sa),  and  C  is  a 
constant  empirically  estimated  to  be  C  •*  3.05  ma/ss.  The  wind  is  assumed  to  be  blow¬ 
ing  in  a  direction  a  with  respect  to  the  x-axls,  and  hence  the  main  ocean  waves  move 
along  this  direction.  The  spectrum  above  is  assumed  non-zero  by  Kinsman  only  in  the 


*Isotropic  means  that  ocean  wave  direction  is  not  recorded  and/or  present.  Some 
directionality  would  be  expected,  with  dominant  ocean  waves  moving  in  the  wind 
direction.  Unfortunately,  too  few  measurements  have  been  made  to  make  a  strong  case 
for  any  given  directional  dependence. 
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half-space  into  which  the  wind  is  blowing;  this  is  reasonable,  for  one  would  expect 
nearly  all  of  the  waves  to  be  driven  in  the  wind  direction.  The  naan-square  height 
of  the  ocean  waves  from  the  spectrum  is  O*  ■  j  C(^/2)a/a  *  (U/2g)#. 

We  now  employ  (44)  in  Equation  (41)  and  evaluate  the  integral  numerically 
for  the  two  extremes  in  wind  (and  wave)  direction:  a  -  0,  representing  radio 
propagation  along  the  dominant  water-wave  motion,  and  a  -  tt/2  representing  propagation 
across  the  dominant  water  wave  direction.  We  call  these  two  cases  the  upwind-downwind 
and  crosswind  directions  respectively.  In  the  integration,  we  divide  (44)  by  two  and 
assuae  the  spectrum  exists  over  all  space  instead  of  over  only  the  forward  wind  half. 
This  is  necessary  because  at  a  given  instant,  the  sea  profile  will  appear  "frozen" 
to  a  radio  wave,  and  it  will  not  be  possible  to  tell  from  this  profile  whether  the 
waves  are  moving  forward  or  backward. 

Figures  4  and  S  show  the  curves  of  the  resistive  (solid  lines)  and  reactive 
(dashed  lines)  portions  of  the  normalized  effective  surface  impedance,  4.  The  sea 
water  is  assumed  to  have  a  conductivity  of  4mhos/meter,  which  is  ths  average  between 
the  winter  and  summer  values  in  the  Atlantic.  The  curves  indicate  that  the  most 
severe  changes  in  surface  impedance  occur  for  a  radio  wave  propagating  in  the  direction 
of  the  dominant  ocean  waves  (l.e.,  across  the  corrugations),  as  one  would  expect. 

These  curves  are  really  only  valid  out  to  frequencies  and  wind  speeds  such  that 
k *0*  «  1,  since  thia  is  one  of  the  assumptions  of  the  analysis.  For  example,  the 
curves  could  be  used  up  to  a  frequency  of  about  30  MHz  with  a  wind  speed  of  23  knots. 

3.  The  Phillips  Spectrum 

Ocean-wave  data  gathered  in  recent  years  have  led  many  oceanographers  to 

search  for  an  alternative  to  the  Neumann-Pierson  spectrum.  Phillips  ^  and  Munk^® 

have  preaented  convincing  evidence  to  show  that  the  spatial  spectrum  should 

behave  as  (pa  +  qa)"a,  instead  of  the  (pa  +  qa)"^*  dependence  of  Equation  (44)  for 

large  pa  +  qa.  More  important,  however,  measurements  show  that  the  lower-end  cutoff 

is  much  more  pronounced  than  the  artificial  exponential  factor  attached  to  the 

49 

Neumann-Pierson  model  in  (44).  Thirdly,  slope  measurements  by  Cox  and  Munk  show 

that  the  spectrun  is  much  closer  to  being  isotropic  than  the  cosine-squared  dlrectlon- 

47 

ality  of  (44).  Hence,  they  suggest  the  following  spectrum  : 

*k da3  •  3/2  c  k*(n/2)3/'3  (U/2g)s 
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Figure  5.  Effective  Surface  Impedance,  Z,  vs  Wind  Speed  and  Frequency. 
Neumann- Pierson  Ocean-Wave  Spectrum  and  Propagation  in 
Crosswind  Direction. 
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"<*•«>  "  n(p^  »BqB)M  (45) 

where  B  “■  0.005.  The  spectrum  la  identically  zero  for  /p3  +  q3  <  g/U3  end  also  in 
the  half-space  from  which  the  wind  is  coming.  The  cutoff  exhibited  by  (45)  is  much 
sharper  than  that  of  (44). 

We  use  (45)  in  a  numerical  evaluation  of  (41),  and  the  results  are  shown 
in  Figure  6.  As  before,  we  divide  (45)  by  two  since  the  "frozen"  surface  profile 
must  produce  a  symmetric  spectrum  about  the  origin.  For  (45),  the  mean-square  height 
is  o3  *  j  B  UVg3.  Since  the  curves  are  meaningful  only  for  kfo3  «  1,  they  can  be 
used  for  frequencies  up  to  about  40  MHz  with  a  wind  speed  of  25  knots.  Again,  the 
conductivity  of  the  sea  is  taken  as  4  mhos/meter. 

In  comparison,  the  upwind-downwind  impedances  from  the  Neumann-Pierson  spec¬ 
trum  are  the  largest.  The  impedances  from  the  Phillips  spectrum  are  reasonably  close 
to  those  for  the  upwind-downwind  direction  with  the  Neumann-Pierson  spectrum.  Those 
for  the  crosswind  direction  of  the  Neumann-Pierson  spectrum  are  the  smallest,  as 
expected.  Until  experimental  evidence  is  obtained,  we  believe  that  the  Phillip's 
impedances  are  perhaps  the  best  estimate  of  the  three. 

F.  Basic  Transmission  Loss  Calculations 
34 

In  the  spirit  of  Wait's  analyses  ,  we  can  use  the  effective  (or  equivalent) 

impedance.  A,  in  the  integrals  for  ground-wave  propagation.  To  facilitate  the 

36 

calculation.  Berry  and  Chrisman  of  ESSA  have  prepared  a  computer  program  which 
solves  this  integral  for  the  field  strength  above  a  spherical  earth.  The  program 
selects  one  of  three  techniques  for  solving  the  integral,  depending  on  whether  the 
observation  point  is  clearly  visible  to  the  source,  In  the  penumbra,  or  in  the  deep 
shadow.  In  the  visible  region,  a  saddle-point  evaluation  is  employed.  In  the 
penumbra,  the  computer  does  a  complex  numerical  Integration  around  a  contour  enclosing 
the  poles  of  the  integrand.  In  the  deep  shadow  region,  the  residue  series  is  used  ;  it 
involves  Fock's  forms  of  the  Airy  functions  and  their  roots  rather  than  the  Hankel 
functions  of  one-third  order,  as  used  by  Bremmer  and  van  der  Pol.  The  computer  output 
has  been  checked  against  Norton's  curves  and  also  against  the  known  solutions  in  the 
visible  region,  and  the  results  agree. 
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We  employ  this  computer  program  to  account  for  sea-state  effects  with  the 
following  minor  modifications:  (i)  we  alter  the  input  format  so  that  surface 
impedance,  A,  can  be  used  rather  than  conductivity  and  dielectric  constant,  and 
(ii)  we  add  another  variable  to  the  output,  basic  transmission  loss,  L^. 

Basic  transmission  loss  is  a  concept  widely  publicized  by  Norton^’ ^  in 
the  1950's.  Formally,  it  is  defined  as 


l*jj  *  10  Logio 


(46) 


where  is  the  power  transmitted  by  an  isotropic  radiator  and  is  the  power 
received  by  an  isotropic  radiator.  In  a  simple  communication  problem,  one  must  merely 
subtract  out  the  free-space  antenna  gains  (in  dB)  in  order  to  determine  the  overall 
power  loss.  For  example,  the  basic  transmission  loss  between  two  points  in  free 
space  separated  by  distance  d  is  10  Logi0  (4nd/ X)3 ;  if  the  same  two  points  are  located 
above  a  flat,  perfectly  conducting  ground  plane,  the  4  is  replaced  by  2  in  the 
parentheses.  The  presentation  of  ground-wave  attenuation  in  the  form  of  basic  trans¬ 
mission  loss  appears  to  be  more  readily  interpretable  and  useful  than  many  of  the 
earlier  ground-wave  normalizations  (e.g.,  unity  electric  dipole,  1  kW  transmitted 
power,  field  strength  one  mile  from  the  transmitter,  etc.).  Examples  illustrating 
the  application  of  these  basic  transmission  loss  curves  will  be  provided  in  the  next 
section. 


For  reference,  we  show  curves  of  basic  transmission  loss  just  above  a  per¬ 
fectly  smooth,  perfectly  conducting  planar  surface  in  Figure  7,  as  a  function  of  range. 
The  values  for  loss  shown  are  essentially  6  dB  lower  than  the  loss  through  free  space 
because  the  vertically  polarized  field  radiated  by  an  isotropic  source  on  the  surface 
is  twice  as  large  with  the  plane  present.  Figure  8  shows  the  basic  transmission  loss 
over  a  perfectly  smooth  spherical  sea  vs  range.  These  curves  were  computed  from  the 
ESSA  ground-wave  program,  using  -  80,  a  -  4  mhos/meter,  and  effective  earth  radius 
“  4/3  actual  earth  radius.  The  former  conductivity  la  typical  of  sea  water  in  the 
Atlantic  ocean.  The  curves  are  calculated  for  frequencies  spanning  the  HF  and  lower 
VHF  regions.  At  short  ranges,  the  curves  of  Figure  8  coincide  with  those  of  Figure  7, 
as  they  should. 


BATTELLE  MEMORIAL  INSTITUTE  -  COLUMBUS  LABORATORIES 


?''S/  WtyLAFlfcnpy 


Figure  7.  Basic  Transmission  Loss  Across  a  Flat,  Perfectly  Conducting,  Perfectly 
Smooth  Ground  Plane  vs  Range  and  Frequency.  Transmitter  and  Receiver 
Antennea  are  at  the  Surface. 
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Figure  8.  Basic  Transmission  Loss  Across  the  Ocean  Between 
Points  at  the  Surface  of  Smooth  Spherical  Earth. 
Conductivity  13  4  mhos/meter  and  an  Effective 
Earth  Radius  Factor  of  4/3  is  Assumed. 
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To  clearly  shew  chs  effects  of  see  state  and  separate  them  from  the  normal 
attenuation  experienced  by  a  ground-wave  above  a  smooth  spherical  sea,  we  compute 
the  difference  between  for  the  smooth  sea,  as  shown  In  Figure  8  and  for  the  rough 
sea.  The  latter  1^  for  the  rough  sea  Is  computed  from  the  ESSA  program  using  the 
values  of  A  (surface  Impedance)  obtained  earlier,  accounting  for  sea  state.  The 
difference  in  transmission  losses  between  the  smooth  and  rough  sea  Is  termed  Ls#, 
and  shown  In  Figures  9  through  32.  Again,  the  conductivity  of  the  sea  water  is 
taken  as  4  mhos/meter  and  the  effective  earth  radius  factor  (accounting  for  refractlv- 
ity)  is  4/3.  The  source  and  receiver  are  both  assumed  to  be  located  at  the  surface. 
Eight  frequencies  are  selected  (3,5,7,10,  15,  20,  30,  and  50  MHz)  which  span  the  HF 
and  lower  VHF  region.  Figures  9  through  16  represent  the  Neumann-Pierson  ocean 
spectrum  for  propagation  in  the  upwind-downwind  direction;  Figures  17  through  24 
use  the  same  spectrum,  but  for  propagation  in  the  crosswind  direction.  Finally  Figures 
25  through  32  present  results  for  the  isotropic  version  of  the  Riillips  ocean  spectrum. 

Figures  33  through  42  show  a  profile  of  basic  transmission  loss  to  various 
points  on  and  above  the  earth's  surface.  The  source  is  assumed  to  be  located  at  the 
surface.  The  first  number  at  each  grid  point  represents  the  transmission  loss  when 
Che  sea  la  perfectly  smooth,  while  the  second  number  represents  the  transmission  loss 
for  Sea  State  5  (25  knot  wind)  using  the  Phillips  spectrum. 

V.  EXAMPLES  OF  THE  USE  OF  THE  GRAPHS 

A.  Surface-to-Surfaca  Conmunicatlon  Problem 

A  shore  station  with  an  array  antenna  is  to  conmunlcate  with  a  ship  300  km 
away  on  a  frequency  of  10  MHz.  The  shore-based  array  has  an  equivalent  free-space* 
gain  of  G£  *  17  dB  and  the  shipboard  antenna  has  a  free-space*  gain  of  *  6  dB 
when  the  two  beams  are  optimally  aligned.  The  shore  station  transmits  an  average 
power  of  1  kW.  We  are  to  find  the  total  power  at  the  receiver  terminals  before 
amplification  for  (a),  a  smooth  sea,  and  (b),  a  sea  generated  by  a  25-knot  wind. 


*This  means  that  the  gains  are  measured  or  calculated  as  though  the  antenna  were 
isolated  in  free  specs  rather  than  being  located  over  a  highly  conducting  ground. 
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Figure  9.  Added  Transmission  Loss  Due  Co  Sea  SCate  at  3  MHz.  Antennas  arc  Located 
Just  Above  Surface.  Neumann-Pierson  Ocean-Wave  Spectrum  with  Prop¬ 
agation  in  Upwind-Downwind  Direction. 


Figure  10.  Added  Transmission  Loss  Due  to  Sea  State  at  5  KHz.  Antennas  are  Located 
Just  Above  Surface.  Neumann-Pier son  Ocean-Wave  Spectrum  with  Prop¬ 
agation  in  Upwind -Downwind  Direction. 
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Figure  11.  Added  Transmission  Loss  Due  to  Sea  State  at  7  MHz.  Antennas  are  Located 
Just  Above  Surface.  Neumann-Plerson  Ocean-Wave  Spectrum  vlth  Prop¬ 
agation  in  Upwind-Downwind  Direction. 
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Figure  12.  Added  Transfusion  Lose  Due  to  Sea  State  at  10  MHr.  Antennaa  are  Located 
Just  Above  Surface.  Neuaann-Plerson  Ocean-Wave  Spectrua  with  Prop- 
agation  in  Upwind-Downwind  Direction. 
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Figure  13.  Added  Transmission  Loss  Due  to  Sea  State  at  15  MHz.  Antennas  are  Located 
Just  Above  Surface.  Neumann-Pierson  Ocean-Wave  Spectrum  with  Prop¬ 
agation  in  Upwind 'Downwind  Direction. 
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Figure  14.  Added  Transmission  Lost  Due  to  Sea  State  at  20  HHs.  Antennas  are  Located 
Just  Above  Surface.  Neumann-Pierson  Ocean-Wave  Spectrum  with  Prop¬ 
agation  in  Upwind -Downwind  Direction. 
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FIGURE  15.  Added  Transmission  Loss  Due  to  Sea  State  at  30  MHz.  Antennas  are  Located 
Just  Above  Surface.  Veumann-Pierson  Ocean-Wave  Spectrum  with  Prop¬ 
agation  in  Upwind -Downwind  Direction. 


FIGURE  16.  Added  Tranamlielon  Lose  Due  Co  See  State  at  SO  KHz.  Antennae  ate  Located 
Juat  Above  Surface.  Neur  jiui-Pierson  Ocean-Wave  Spectrum  with  Prop¬ 
agation  in  Upwind-Downwind  Direction. 
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Figure  17.  Added  Transmission  Loss  Due  Co  Sea  State  at  3  MHz.  Antennas  are  Located  Just  Above  Surface 
Neumann- Pierson  Ocean-Wave  Spectrum  with  Propagation  in  Crosswind  Direction. 
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Figure  18.  Added  TransnUaion  Lore  Due  to  Sea  State  at  5  MHz.  Antenna*  are  Located  Juat  Above  Surface 
Neunano-Plerson  Ocean-Wave  Speccrun  with  Propagation  in  Croaswlnd  Direction. 
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Added  Transmission  Loss  Due  to  Sea  State  at  7  MHz.  Antennas  are  Located  Just  Above  Surface 
Neumann-Pf erson  Ocean-Wave  Spectrum  with  Propagation  in  Crosswind  Direction. 
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Figure  2U.  Added  Transmission  Loss  Due  co  Sea  State  at  10  MHz.  Anteonaa  are  Located  Juat  Above  Surface 
Neumann- Pier eon  Ocean-Wave  Spectrum  with  Propagation  In  Croaevlnd  Direction. 


Figure  21.  Added  Transmission  Loss  Due  to  Sea  State  at  15  MHz.  Antennas  are  Located  Just  Above  Surface 
Neumann- Pierson  Ocean-Wave  Spectrum  with  Propagation  in  Crosswind  Direction. 
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Figure  22.  Added  Transfusion  Lose  Due  to  Sea  State  at  20  MHz.  Antennae  are  Located  Juat  Above  Surface 
Neunann- Pier son  Ocean-Wave  Spectrins  with  Propagation  in  Crosawind  Direction. 
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FiGURK  23.  Added  Transmission  Loss  Due  to  Sea  State  at  30  MHr.  Antennas  are  Located  Just  Above  Surface. 
Neuaann-Pierson  Ocean-Wave  Spectrum  with  Propagation  in  Crosswind  Direction. 
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FIGURE  24,  Added  Transmission  Lose  Due  to  See  Stete  et  50  MHz .  Antennas  ere  located  Just  Above  Surface 
Neumann-Pierson  Ocean-Wave  Spectrum  with  Propagation  in  Crossvind  Direction. 


Figure  25.  Added  Transmission  Loss  Due  Co  Sea  SCate  at  3  MHz.  Antennas  are  Located  Just  Above  Surface 
Phillips  Isotropic  Ocean-Wave  Spectrum. 
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Figure  26.  Added  Transmission  Loss  Due  Co  See  State  at  5  KHz.  Antennas  are  Located  Jual  Above  Surface 
Phillips  Isotropic  Ocean-Uave  Spectrum, 
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Figure  27.  Added  Transmission  Loss  Due  Co  Sea  State  at  7  KHz.  Antennas  are  Located  Just  Above  Surface 
Phillips  Isotropic  Ocean-Wave  Spectrum. 
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Figure  28.  Added  Transmission  Los*  Due  to  See  State  at  10  HHi.  Antennae  are  Located  Juat  Above  Surface 
Phillips  Isotropic  Ocean-Wave  Spectrum. 


Added  Transmission  Loss  Due  to  Sea  State  at  15  MHz.  Antennas  are  Located  Just  Above  Surface 
Phillips  Isotropic  Ocean-Wave  Spectrum. 
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Figure  30.  Added  Transmission  Lost  Due  Co  See  State  at  20  MHz,  Antennas  are  Located  Juat  Above  Surface 
Phtlllpa  Isotropic  Ocean-Wave  Spectrum. 
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Figure  31.  Added  Transmission  Loss  Due  to  Sea  State  at  30  MHz.  Antennas  are  Located  Just  Above  Surface 
Phillips  Isotropic  Ocean-Wave  Spectrum. 
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Figure  32.  Added  Tranrnlsalon  Loe*  Due  to  Sea  State  at  50  MHi.  Antenna*  ar*  Located  Juat  Above  Surface 
Phillip*  Isotropic  Ocean-Wave  Spectrun. 
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figure  Jfc.  gaelc  Transmission  Lose  to  Points  «t  Various  Height*  and  Mug  •»  Above  th*  Oct  mb  at  }  MHa. 
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The  total  loss  la  given  as 

pt 

10  Logi0  p5,  “  S  “  Gt  ‘  Gr  ‘  <47> 

From  Figure  8,  we  find  that  at  300  km,  «  121.8  dB  for  a  smooth  sea;  uaing  Figure  28 

we  note  that  this  increases  to  132.7  dB  for  a  sea  driven  by  a  25-knot  wind.  Hence, 

L.  -  G„  -  G_  *  98.8  dB  and  109.7  dB  for  these  two  sea  conditions.  Thus,  for  P  •  1  kW 
d  c  r  t 

(average),  1.3  X  10"7  Watts  (average)  when  the  sea  is  calm,  but  drops  to  P^ 

1.1  X  10“*  Watts  (average)  in  rough  seas  produced  by  a  25-knot  wind. 

B.  Surface-to-Air  Communication  Problem 

Here  let  ua  consider  a  ship-to-air  communication  application.  The  frequency 
is  20  MHz .  We  wish  to  determine  the  signal  power  received  by  an  aircraft  at  500  meters 
altitude  (about  1500  feet)  and  270  km  away  from  the  ship  during  calm  sea  conditions' 
for  500  watts  transmitted  power.  Furthermore,  we  wish  to  estimate  how  much  higher 
the  plane  must  fly  *t  that  range  to  receive  the  same  signal  power  when  the  sea  is 
aroused  by  a  25  knot  wind.  Possible  ionospheric  reflections  are  not  to  be  included 
in  the  signal  estimates.  The  aircraft  free-space  antenna  gain  is  about  6  dB. 

For  the  smooth-sea  problem,  we  employ  Figure  40  to  find  the  basic  trans¬ 
mission  between  two  isotropic  sources,  one  located  on  the  surface  and  the  other  at 
500  meters  altitude  and  270  km  range.  This  is  given  by  the  first  number  at  the 
appropriate  grid  point  on  the  chart,  i.e.,  -  145.7  dB.  Then  using  Equation  (47), 

we  find  that  the  received  power  will  be  Pf  **  1.35  X  10“11  watts. 

At  that  same  altitude,  we  see  from  the  second  number  at  the  same  grid  point 
that  the  transmission  loss  increases  to  153.9  dB  when  the  sea  is  fully  aroused  by  a 
25  knot  wind.  This  8.2  dB  drop  in  received  signal,  however,  can  be  recovered  by 
climbing  in  altitude  to  the  next  grid  point  (i.e.,  1000  meters  or  3000  feet),  at  which 
the  loss  drops  back  to  145.1  dB,  near  the  original  145.7  dB  figure  for  the  smooth 
sea. 
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VI  ■  SUMMARY 

This  report  hes  had  one  principal  purpose  throughout:  to  derive  some 

estimates  of  the  effects  of  sea  state  on  ground^rave  propagation  at  HP  and  VHP. 
f  % 

At  the  outset  of  the  analysis  it  was  established  that  the  only  truly  applicable 
boundary  condition  is-  the'  Laontovich,  or  impedance  boundary  condition.  This 
boundary  condition,  when  applied  to  the  sea,  correctly  accounts  for  the  slight  losses 
associated  with  the  finite  water  conductivity.  With  this  boundary  condition,  the 
derivation  followed  the  classical  statistical  boundary  perturbation  approach 
developed  by  Rice.  The  effective  surface  Impedance  accounting  for  the  roughness 
was  derived  in  this  manner  for  s  vertically  polarized  wave  which  appears  to  propagate 
locally  over  the  surface  at  the  pseudo-Brewster  angle.  Since  the  ground  wave 
propagates  locally  at  this  angle,  the  effective  surface  impedance  thus  derived  can 
be  used  in  any  of  the  conventional  treatments  of  radiation  and  propagation  over  a 
spherical  earth  with  an  Impedance  boundary. 

The  effective  surface  Impedance  obtained  here  is  seen  to  consist  of  two 
terms:  one  which  is  She  impedance  of  a  otherwise  smooth  planar  surface  of  aea  water, 
and  the  second  which  is  due  to  the  roughness  height  spectrum  of  the  ocean  waves. 

This  latter  term  preserves  the  conductivity  dependence,  and  the  error  involved  in 
assuming  a  perfectly  conducting  surface  can  thus  be  clearly  seen.  The  increase  in 
the  surface  impedance  due  to  the  roughness  is  due  largely  to  two  phenomena:  the 
resistance  increase  is  caused  by  removal  of  energy  by  Bragg  scatter  from  the  ground 
wave.  The  inductive  reactance  increase  is  due  to  the  shorter  roughness  wavelengths 
which  do  not  scatter,  but  produce  evanescent  modes  near  the  surface.  The  latter 
phenomenon  is  identical  to  that  behind  the  use  of  short-period  corrugations  to  "trap" 
energy  near  guiding,  slow-wave  structures. 

The  analysis  presented  ha a  emphasized  the  electromsgnetlc  derivation. 

No  attempt  was  made  to  analyze  or  improve  the  oceanographic  models  for  the  sea  height 
spectrum.  In  order  to  obtain  estimates  of  sea  state  effects,  two  coesaonly  seen 
ocean-wave  spectra  models  were  used:  the  Neumann- Pier son  and  the  Phillips  spectrum. 

The  effective  surface  impedance  based  on  these  two  models  was  numerically  calculated 
for  various  frequencies  and  wind  conditions.  Finally,  we  used  these  surface  Impedances 
in  a  standard  program  to  compute  the  basic  transmission  loss  across  the  sea  vs  range 
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for  various  frequencies  and  winds.  Of  the  two  ocean-wave  models,  it  is  believed  that 
the  Riillips  spectrum  is  somewhat  more  Justified  than  the  Neuaann-Pierson  directional 
model. 

To  illustrate  how  these  curves  can  be  applied  to  connunication  problems  we 
traced  through  some  sample  problems.  The  curves  can  also  be  applied  to  radar 
applications  in  a  straightforward  manner. 

VII.  CONCLUSIONS 

Wind  speed  and  the  resulting  ocean  roughness  will  definitely  affect  the 
propagation  of  a  ground  wave  in  the  HF/VHF  regions.  This  roughness  effect  is 
generally  negligible  below  a  frequency  of  3MHz.  At  15  MHz  and  a  range  of  100  miles, 
the  received  signal  over  the  Atlantic  is  shewn  to  vary  as  much  aa  15  dB  due  to  sea 
state.  Above  about  50  MHz  and  for  wind  speeds  greater  than  25  MHz,  the  techniques 
derived  here  are  no  longer  valid  because  of  the  assumptions  inherent  in  the 
analysis. 

Any  deviations  of  measured  propagation  loss  from  that  predicted  by  our 
curves  are  not  believed  due  to  any  shortcomings  of  the  electromagnetic  analysis. 

They  are  due  rather  to  the  empirical  ocean-wave  spectra  models  used  here.  Reliable 
directional  spectra,  as  required  by  our  equations,  have  not  been  measured  in 
sufficient  detail  by  oceanographers  to  permit  the  establishment  of  an  empirical 
model.  In  addition,  the  only  isotropic  spectra  models  available  are  based  on  wind- 
driven  seas  and  do  not  Include  the  effects  of  swell  generated  by  storms  in  other  areas. 
There  is  at  present  no  reasonable  estimate  of  the  error  resulting  from  this  neglect  of 
swell. 

Measurements  of  ground-wave  losses  over  the  sea  are  in  progress  and  results 
should  be  available  for  comparison  in  the  near  future.  As  these  measurements  are 
being  made,  Raytheon  is  also  recording  the  sea  conditions  for  correlation  with  the 
signal  strength.  Along  with  their  measurements  of  radar  sea  clutter,  this  data  could, 
as  a  by-product,  provide  welcome  information  about  the  ocean-wave  spectrum. 
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APPENDIX 


Ground-Wave  Bibliography 

The  aubjecc  of  ground-wave  propagation  acroaa  the  earth  is  one  of  the  most 
widely  treated  technical  subjects  of  this  century.  We  attempt  to  list  hare  moat  of 
tha  principal  open-literature  publication  relating  to  this  subject.  We  apologias 
for  omissions,  of  which  there  are  certainly  many.  Tha  list  hare  has  invariably 
failed  to  Include  many  reports,  even  of  an  unclassified  nature,  as  well  ss  many 
foreign  articles  on  the  subject. 

The  bibliography  does  not  Include  treatments  of  tha  effects  of  the  ionosphere 
or  atmosphere.  Nor  is  ground  scatter  (cluttar)  considered. 

So  as  to  show  tha  historical  evolution  of  tha  subject,  wa  break  the  biblio¬ 
graphy  into  decades.  All  of  the  publications  generated  within  a  given  decade  are 
then  arranged  alphabetically  by  author. 

Before  1911 

Epstein,  P. ,  "The  Propagation  of  Waves  in  Wireless  Telegraphy,  Taking  Into  Consider- 
atlon  the  Nature  of  the  Ground",  ^ahrbach  d.  Drahtloaen  Telcgraphie.  Vol  4,  No  2, 
pp  176-187  (December  1910). 

Hack,  F.,  "The  Propagation  of  Electromagnetic  Waves  Over  a  Plana  Conductor",  Ann. 

Rws .  Lpz . .  Vol  27,  p  43  (1908). 

Koepsel,  A.,  "Role  of  the  Earth  in  Wireless  Telegraphy",  Dingier’ s  Polvtechn.  Journal. 
Vol  318,  pp  385-388  (June  20,  1903). 

Nlckolaon,  J,  W.,  "On  the  Banding  of  Electric  Waves  Round  the  Earth",  Part  I,  Phil. 
Mag..  Vol  18,  p  757  (1910). 

Nlckolaon,  J.  W.,  "On  the  Bending  of  Electric  Waves  Round  the  Earth",  Part  11,  Phil. 
Mag..  Vol  20,  p  157  (1910). 

Soumerfeld,  A.,  "The  Propagation  of  Waves  in  Wireless  Telegraphy",  Ann.  Physlk.  Vol  28, 
p  665  (1909). 

Sonmerfeld,  A.,  "The  Propagation  in  Wireless  Telegraphy",  Jahrbach  d.  Drahtloaen 
Telcgraphie.  Vol  4,  p  157  (1910). 

Sotmnarfeld,  A.,  "Uber  die  Fortpfanzung  electrodynamiacher  Wellen  langs  elnes 
Drahtea",  Phyajk  u.  Chem..  Vol  67,  p  233  (1899). 
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Zenneck,  J.,  "Uber  die  Fortpf lenzung  ebener  electromagnetisher  Wallen  langs  einer 
ebenen  Leiterflache  und  ihrer  Beziehung  zur  drahtlosen  Telegraphie",  Ann.  Phvalk. 

Vol  23,  pp  846-866  (1907). 
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